Laser stabilization using spectral hole burning 
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We have frequency stabilized a Coherent CR699-21 dye laser to a transient spectral hole on the 
606 nm transition in Pr^^:Y2Si05 . A frequency stability of 1 kHz has been obtained on the 10 \xs 
timescale together with a long-term frequency drift below 1 kHz/s. RF magnetic fields are used to 
repopulate the hyperfine levels allowing us to control the dynamics of the spectral hole. A detailed 
theory of the atomic response to laser frequency errors has been developed which allows us to design 
and optimize the laser stabilization feedback loop, and specifically we give a stability criterion that 
must be fulfilled in order to obtain very low drift rates. The laser stability is sufficient for performing 
quantum gate experiments in Pr^+:Y2Si05 . 



I. INTRODUCTION 



Quantum information processing and the ability to ex- 
perimentally manipulate quantum systems has evolved 
significantly over the past decade, and the search for 
new and potentially scalable qubit systems is ongoing. 
One of these schemes employs rare-earth-ion-doped crys- 
tals [U, 0, H, 0, where, among others, the materials 
Pr^+:Y2Si05 and Eu^+:Y2Si05 are good candidates with 
very good coherence properties [1, B B @- Qubits are 
represented as hyperfine levels |0) and |1), and they are 
manipulated using an optical transition to an excited 
level |e). This optical manipulation requires very sta- 
ble lasers in order to obtain high-fidelity gate operations. 
In order not to be limited by laser performance, the laser 
must remain phase coherent on a timescale comparable 
to the optical coherence time, T2, of the atomic material. 
Our experimental focus is on Pr^+:Y2Si05 which has an 
optical transition at 605.977 nm with coherence time as 
long as T2 — 152 \\s [7], which is accessible with dye 
lasers only. The narrow optical transitions are hidden in- 
side a much broader inhomogeneous profile, but they can 
be isolated by optical pumping methods [l3, [HI , taking 
advantage of the very long hyperfine level lifetime (being 
ftilOO seconds in Pr^+:Y2Si05 and much longer in a mod- 
erate magnetic field). These initialization techniques em- 
ploy several cycles of population transfer with subsequent 
relaxation and may last for times of the order of several 
hundreds of milliseconds. This also places some long- 
term stability requirements on the laser performance. 

Frequency stabilization of lasers is a science in itself. 
Previously, dye lasers have been stabilized using optical 
cavities E 111 Q [H, H 113, El, H H HI, and re- 
cently the method using spectral hole burning has been 
demonstrated with semi-conductor laser s [22l . |23| . [2I . |25| . 
[11,113, Hi and with a Tirsapphire laser [iJ^We demon- 
strate for the first time a dye laser stabilized to a spectral 
hole. The method employing spectral holes has a number 
of advantages over optical cavities. Firstly, the sensitiv- 
ity to vibrations requires that the atomic medium moves 
much less than an optical wavelength in an optical coher- 
ence time. However, for a high-finesse cavity this sensi- 
tivity is essentially multiplied by the average number of 



round trips made by a photon in the cavity, which can be 
several orders of magnitude. In addition, using the same 
hole-burning crystal for laser stabilization and for further 
experiments may reduce the sensitivity to vibrations sig- 
nificantly. Secondly, to make a phase-coherent laser on 
the timescale of T2 requires a physical reference system 
with a phase memory of comparable duration. Hence, 
when performing experiments on e.g. Pr^+:Y2Si05 , with 
T2 = 152 |is, the material itself will automatically be suf- 
ficient for laser stabilization. A similar round trip time in 
a cavity is possible, but technically very challenging. A 
cavity also requires mode matching, and optical isolators 
are required to prevent the directly back-reflected field 
from entering the laser cavity. The main drawbacks of 
spectral holes include the requirement of cryogenic cool- 
ing, and for transient hole systems the long-term stability 
presents a challenge since the spectral hole position may 
change over time. 



In the following sections we present a detailed theory 
of the dynamics of spectral holes and on the response 
of the atomic medium to errors in the laser frequency. 
This helped us to design the feedback hardware and in 
realizing the laser stabilization experimentally. In doing 
so, we also introduced an "RF-eraser", which consists 
of RF magnetic fields controlling the hyperfine level life- 
times. In this manner, we can adjust and optimize the 
spectral hole dynamics for best laser performance. We 
also demonstrate experimental support for the theoret- 
ical calculations with emphasis on the problem of laser 
drift. The RF magnetic fields are very useful in this con- 
text. Apart from the theoretical understanding, it is also 
very important to build fast, low-noise, and low-offset 
feedback electronics. We have designed analog electronic 
circuits using the best components available on the mar- 
ket for the fastest parts of the regulation system. The 
most important design considerations are given in this 
paper, while the technical details are available elsewhere 

mm- 
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II. THEORETICAL DESCRIPTION OF ION 
LOCKING 

The entire theoretical understanding of laser stabiliza- 
tion to spectral holes in inhomogeneous broadened tran- 
sitions is one of the key results of this paper and is in- 
troduced in this section. Section III Al presents simple 
arguments to explain the basic idea and to compare the 
procedure of locking to spectral holes with the method of 
locking to optical cavities. We then proceed to Sees. Ill Bl 
and III CI where we introduce the Maxwell-Bloch equa- 
tions, written in a form convenient for inhomogeneously 
broadened transitions and then calculate the propagation 
effects of laser light through the atomic material in the 
presence of spectral holes. With these results we are able 
to calculate the error signals suitable for the laser sta- 
bilization scheme using the Pound-Drever-Hall method 
[12] in Sec. Ill Dl This section provides a general under- 
standing of the response of atoms to errors in the laser 
frequency, which will put us in a position to design the 
electronics hardware for the stabilization feedback loop. 
In Sec. lIIEl we analyze the possibility of linear drift errors 
and their prevention, and finally in Sec. Ill Fl we summa- 
rize our theoretical calculations. 

The intention of this paper is to enable the reader to 
understand the details of laser stabilization to spectral 
holes with special emphasis on the feedback system de- 
sign. Experimentally, we focus on a particular atomic 
system consisting of ions in a Pr^+:Y2Si05 crystal. In- 
stead of developing a general theory we make a few ap- 
proximations suitable for this particular system. How- 
ever, we present the calculations in such detail that it 
should be clear how to extend the theory to other sys- 
tems. 

In order to maintain the physical understanding we will 
restrict ourselves to analytical derivations and make ap- 
proximations, rather than numerical simulations, when 
the calculations become difficult. Our theory is quanti- 
tatively accurate for many practical systems. 

For readers not interested in all the technical details 
we will now summarize the most important theoretical 
results. After the simple discussion of laser stabiliza- 
tion schemes in Sec. Ill Al we proceed to Sec. Ill Bl where 
we state in Eqs. (|3]|ni) the important propagation equa- 
tions for different frequency components of the laser light 
through matter. In Sec. lIICl these are applied to the spe- 
cial case of a spectral hole in a broad inhomogeneous pro- 
file in Eqs. and (pij) which, together with the hole 
shape parameters of Eqs. (|14m6p . form the workhorse 
for all calculations of error signals. In Sec. Ill Dl the error 
signals are calculated; the main result being the trans- 
fer function (|36p connecting errors in laser frequency to 
modulations in the detected laser power. These calcu- 
lations also indicate the optimum parameter choices for 
best performance. To mention one example, the modula- 
tion index m = 0.56 in general optimizes the magnitude 
of the error signal which is different from the optimum re- 
sult m = 1.08 known from cavity locking 32]. However, 



additional constraints on the modulation index may be 
necessary when we take the possibility of laser drift into 
account in the rather technical section Hi El Here we cal- 
culate the correction ([55)1 to the error signal if the laser 
is drifting linearly with time. This specifically leads to 
a criterion (j59p which must be fulfilled in order to pre- 
vent laser drift. A practical version of this criterion is 
illustrated in Fig. [31 



A. Stabilization - an overview 

Before entering into the detailed theoretical description 
of laser locking to spectral holes, we will give a general 
picture of the stabilization mechanisms. We discuss the 
concept of frequency locking and phase locking, and to 
this end we start by briefly discussing how stabilization to 
a cavity works. An extensive and pedagogical description 
can be found in [s^ . 

In stabilization setups with optical cavities the laser 
light is sent into a high-finesse cavity and the back- 
reflected light is detected. The aim is to maintain the 
laser frequency at the cavity resonance at all times, in 
which case the reflected field is ideally zero (all light is 
transmitted through the cavity). Let us assume that the 
cavity mirrors reflect 99.99% of the light intensity. The 
fact that no light is reflected in the ideal case is due to 
the fact that the immediately reflected laser field at the 
incoming mirror interferes destructively with the 0.01% 
leakage of the 10''' times more intense intra-cavity field. 

Let us now assume that the incident light suddenly 
changes phase on a timescale much shorter than the cav- 
ity build-up time, r = ^ (here T — 0.01% is the mirror 
transmission, c is the speed of light, and L the cavity 
length. For L = 10 cm we have r « 3 The imme- 
diately reflected beam then changes its phase and does 
not interfere destructively with the cavity leakage field, 
which retains its original phase for a time of the order 
of r. The intensity is then non-zero and an error can be 
detected. The immediate phase difference (including the 
sign) of the two fields can be measured using the Pound- 
Drever-Hall method [l2| • If the phase change persists for 
a time longer than r the reflected field returns to zero. 

We also consider the case where the frequency is 
changed from the cavity resonance. On a timescale longer 
than T this field will build up inside the cavity, but since 
the resonance condition is not met, the cavity length is 
not equal to an exact integer multiple of A/2, where A 
is the wave length of the light. Hence, each round trip 
in the cavity is slightly advanced or delayed in phase, 
and the total contribution will be phase shifted, with a 
sign depending on the sign of the frequency error. When 
this field leaks out of the cavity and interferes with the 
immediately reflected field, we can again, in a similar 
fashion to above, measure the frequency error (including 
the sign). Note, that the non-zero reflection persists over 
time and is not limited to a duration r. 

The above two examples illustrate the fact that on long 
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timescales the detection and correction of frequency er- 
rors is possible, while phase errors can only be detected 
for a limited time given by the phase memory time t. 
We shall now argue that this picture also holds for lock- 
ing to spectral holes. In this scenario the laser light 
is incident on an inhomogeneously broadened transition 
and the laser will burn a spectral hole in the absorbing 
background. To some extent, this spectral hole can be 
thought of as the inverse of the pedagogically simpler 
case with no background absorption and the presence of 
atoms at a single frequency only. We adopt this picture 
below. 



It is well known that when scanning a laser across an 
atomic resonance the absorption will change (proportion- 

ally to -J- / + 1\ ) where Fh is the full width half max- 
imum (FWHM) line- width and A the detuning). At the 
same time, the index of refraction will change propor- 
tionally to A^^/(^ -I- A^) which, for small detunings 
A, causes a phase shift of the light proportional to the 
detuning. This effect can be detected by methods simi- 
lar to the Pound-Drever-Hall method [I4I and is used for 
measuring frequency errors. 

The above picture describes a steady state which re- 
quires modifications on short time scales. If the phase 
of the incoming light changes abruptly, it takes of the 
order of the optical coherence time T2 = 7?- before the 

i h 

atoms have reached a new steady state (Fh is here mea- 
sured in rad/sec). Meanwhile, the atoms keep radiating 
at the original phase, which allows the phase error to be 
measured but only for a limited time, T2. 

We see that the methods of optical cavities and spec- 
tral holes are very similar. In a common picture we can 
think of both as frequency filters where the phase mem- 
ory is set by the inverse line-width and the filtering gives 
rise to slowly varying fields. The possibility of detect- 
ing phase errors relies on the interference between an 
immediately changing field and the slowly varying field. 
For the optical cavity the immediately changing field is 
the directly reflected beam, while for spectral holes it 
is the remaining part of the incoming light transmitted 
through the atoms. The slowly varying field is the intra- 
cavity field in the case of optical cavities, while it is the 
light emitted coherently by atoms in the case of spectral 
holes. The Pound-Drever-Hall method [12] allows for the 
detection of the above mentioned interference in both 
cases, but there is one important difference. The modu- 
lation sidebands required for this scheme are not stored 
in the optical cavity (they are not resonant) and hence do 
not experience the frequency filtering. For the method 
of spectral holes these sidebands create side holes which 
must be taken into account to understand the complete 
picture. 



B. Two-level atoms and Maxwell-Bloch equations 

For the detailed theoretical calculations we start with 
an ensemble of inhomogeneously broadened two-level 
atoms. We allow laser light to propagate through these 
along the z-direction. With a large beam cross section 
a one-dimensional theory is sufficient, and the Maxwell- 
Bloch equations can be written (see e.g. [H, [H, [sH] ) : 

d V\ 

— (m — iv) = ^(~2~ ^ — iv) — iflw, (1) 

^u; = ^ [n{u + iv) ^n*[u^ iv)] - ^(1 + w), (2) 

Here {u, v, w) is the usual Bloch- vector which depends on 
time i, position z, and detiming A (from a chosen refer- 
ence point). The electric field is described in terms of 
the complex Rabi frequency Q{z,t) = ii£{z,t)/h, where 
£ is the complex electric field and /i is the dipole moment 
along the direction of the field (we consider only a single 
linear polarization mode). Fh — 2/T2 is the FWHM ho- 
mogeneous line- width of the atoms in rad/sec, Ti and T2 
are the life and coherence times of the optical transition, 
respectively. In Eq. ([3]) n is the refractive index of non- 
absorbing background atoms, and (?(A) is a dimension- 
less function describing the inhomogeneous distribution 
of atoms such that g{A) is proportional to the number of 
atoms with transition frequency A. We use an unconven- 
tional but experimentally convenient normalization such 
that if the absorption length measured with a weak laser 
field at frequency Aq amounts to ao, then g{Ao) must 
be unity. The integral over A in Eq. ^ effectively adds 
the contribution of the polarization from all the atoms to 
the electric field J7 at position z and time t. 

Eqs. ((T][3]) are in general difficult to solve analytically. 
However, for our specific needs regarding laser stabiliza- 
tion we will make a number of approximations in the 
following. We start by noting that in Eq. ^ the term 
is only relevant when describing very fast changes, 
on the time scale L/c where L is the length of the sample. 
We therefore neglect this term completely. 



1. Linear regime of Maxwell-Bloch equations 

The next approximation is to consider Eqs. (Il][3]) in the 
linear regime where, i/; « — 1 for all atoms, i.e. the prob- 
ability of being in the excited state is small. In Sec. Ill CI 

we discuss the validity of this approximation. Insert- 
ing (with w — —1) the integral form u{z,t) — iv{z,t) — 

i Jl^e'^^+'^'>^*'*">n{z,t')dt' ofEq. (P) intoEq. dll) and 
expressing the electric field fl in terms of its Fourier com- 
ponents, ^l{z,t) = J_^il{z,uj)e~^'^*duj, it follows that 
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Eq. ([3]) can be written in Fourier space as: 



d_ 

dz 



Lo) 



'27r 



g(A)rfA 

-oo +i(A-cj) 



Uj) 



aR(w) + iQ;i(cj) 



(4) 



where we define a-{i{Lj) and ai(aj) as the real and imag- 
inary absorption lengths, respectively. For a single fre- 
quency component of the field n{z,oj) = A{z,Lu)e~^'^^^-''^^ 
with real amplitude A and phase (f> we have the relation: 



dA{z,Lu) 

dz 
d4>{z, uj) 

dz 



2 



Aiz,Lu), 



(5) 
(6) 



Eqs. ([Uini) will be the workhorse for many calculations in 
the following sections. Our goal is to model the frequency 
variations of the incoming laser field, propagate this field 
through the atomic medium via Eqs. ([5]) and ([6]), and 
finally derive an error signal useful for frequency stabi- 
lization based on the outgoing field. We will reach this 
goal in Sec. Ill Di but before that we introduce a model 
which describes the effect of hole burning in terms of the 
shape function, (?(A). 



2. Example: Inhomogeneous profiles 

Let us give a simple and useful example of the above 
equations. Consider a Lorentzian-shaped inhomogeneous 



profile, g{A) 



I . Using the residue 



theorem we find (for Finh ^ Eh): 



|e) states we add a third reservoir state |r) (see Fig.|4l[b)) 
and write the Bloch equations for these: 



d_ 
di 



{u - iv) = -{-z- + iA){u ~ iv) - iVL{pe - Pg), (8) 



dpe 

dt 
dt 



dpr 

dt 



'-[n{u + iv)-n*{u-iv)]-^^pe. (9) 



— — [Vl{u + iv) ~ Vl*{u — iv) 



-Pa 



-Pr, 



■Pe 



Pa 



-Pr- 



(10) 

(11) 



We assume the optical field, fJ, only couples to the tran- 
sition |g) |e) and hence u and v still refer to this 
transition, and only the factor w = pe — Pg appears in 
the driving term in Eq. ([5]), as was the case in Eq. ([1]). 
We add the possibility of decays from the excited state 
to the reservoir state. The branching ratios from |e) to 
|g) and |e) to |r) are denoted beg and 6or, respectively. 
We also model relaxation between the |g) and |r) levels. 
The timescale for decays from |g) to |r) is Tgr, which in 
general need not be the same as the timescale Tj-g in the 
opposite direction. For the homogeneous line- width, Fh, 
we now have: 



Eh 
2 



1 



1 



1 



(0) 



IT ' 



(12) 



where Tj*'^ is the coherence time of the optical transition 
|g) — > |e) in the absence of ground state relaxation, and 
the term takes the finite lifetime of the state |g) into 
account. 



aR(w) = ao 



r?„, 



(7) 



1. Separation of times cales 



Here we have chosen the frequency zero point to be at the 
center of the profile. We see that the absorption length at 
this frequency is ag, consistent with g(A) — 1 for A = 0. 

We could also consider a Gaussian-shaped profile, 
g(A) = exp(-41n(2)AVFf„i^), with FWHM Finh- In- 
serting this into Eq. Q we find with help from [36j . 



53 ~ 



aow(Z) 



Q{z,uj), where Z — 



2\/ln 2 / 



2 ao^^,^;, ^ — -^^—^{Uj + i-^) 

and w{Z) is the error function for complex arguments. 
Thus q;r = QoRe|w(Z)| and a\ — aoIm{w(Z)}, known 
as the Voigt profile [3^ . In general, the shape of the in- 
homogeneous profile varies depending on the broadening 
mechanism 13 71. 13811 . 



C. Two-level atoms with a reservoir state 

The calculations in the previous sections need to be 
refined in order to describe the effect of spectral hole 
burning. So, in addition to the ground |g) and excited 



Our next step from Eqs. ([SlfTTjl is to use steady-state 
solutions for pe, Pg, and p-c while u and v still are allowed 
to vary in time according to Eq. ([8]). This is a good ap- 
proximation since in our specific case we have naturally 
different timescales for the ground state populations and 
the optical coherence, Tj-g, Tgr 12- Furthermore, when 
the laser is actively stabilized to a line-width narrower 
than Fh, it is a good approximation to assume a zeroth 
order starting point, Vl = fioe"*^"*, where the laser is 
running perfectly at a monochromatic frequency, Aq. If 
the variations from this starting point are small, the pop- 
ulations will always be close to their steady-state values. 



So, with Q. = Jlpe 



-iAnt 



we calculate the steady-state 



solutions of Eqs. (|5]fTT|l. With a little work we obtain 
and expression for the population difference pg — po {pc, 
Pg, and Pr are given in detail in Eq. ([54]) ): 



— Pc = G \ 1 — rfholc - 



+ (A - Ao) 



(13) 
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where dhoic is the relative hole depth and Fhoie is the 
FWHM of the hole. These parameters can be written: 



Y^lT+^Rj^' Thole = Thjl + (1 + i?)^, 



So 
2 

(14) 



where sq is the resonant saturation parameter: 

so = If^opTiTa, (15) 
and for our particular case of Eqs. ((8l fTT|) we have: 



R 



1 + ^ 1 



Ti 



1 + 



Zk 



G 



1 + 



Zk 



(16) 



The saturation parameter so is a measure of the probabil- 
ity of an atom being in the excited state |e) at resonance. 
In steady state at A = Ao we have ^ = (l + f-). 
The parameter G is a measure of the fraction of atoms 
in the ground state |g) in equilibrium in the absence of 
the laser light or when the detuning, A — Aq, is large. 
The parameter i? is a measure of how likely it is for an 
atom to be trapped in the reservoir state. The essence of 
R is in the term ^crTi-g/Ti, which is the ratio of the rate 
bei/Ti from |e) into the reservoir state |r) and the rate 
1/Trg out of the reservoir state. From Eq. (fT4|) it is clear 
that a spectral hole can be deep and broad for different 
reasons: Firstly, if the laser field is strong with a high 
saturation parameter so, although R is small, and sec- 
ondly, if the trapping parameter R is large even a weak 
field with sq ^ 1 is capable of digging a deep, wide hole. 

Now, we wish to employ Eq. ([3]) or (jl]) together with 
Eqs. PfTTj) . We insert the steady-state value of Eq. (fTB]) 
into Eq. ([5]) written in integral form. Since the steady- 
state value is time- independent we may perform the same 
steps as those leading to Eq. We will incorporate 
the value of pg — pc into the g{A) shape function and 
just pretend that we never left the linear approximation, 
Eq. ([4]), of a two-level system. This is done correctly 
when: 



5(A) 



Pg 



^holc 4 



G 



(A-Ao)2 



(17) 



equations (|4]) and ((T7)) . The (/(A) function for a spectral 
hole describes how many atoms actually participate in 
the active two-level transition |g) |e). We included 
the effect of saturation where atoms can also populate 
the excited state |e) (which mathematically also creates 
a hole in pg — pc). However, since we assumed the popula- 
tions Pc, Pg, and pr to be essentially constant in time, we 
have restricted ourselves to solutions where the laser field 
does not deviate much from a perfect field, ft = iloe"'^"* 
(we have linearized the theory around this zeroth order 
solution). Note, that the field can still have fast varia- 
tions in e.g. its phase, as long as the phase excursions 
are not too large. Since both population trapping in 
the reservoir state and the effect of saturation (leading 
to population trapping in the excited state) are incor- 
porated into the single parameter (/(A), we effectively 
model the three-level equations ([8l lTT|) with our initial 
linear two-level system with low saturation, as described 
by Eq. ©. 

Using Eq. (|17p also requires another approximation. 
We note that dhole and Fhoie depend on the resonant sat- 
uration parameter, sq. If the optical depth, agL, of the 
atomic sample is large, the saturation parameter will de- 
pend on z, and the use of a z-independent 5(A) will be 
incorrect. However, if the laser field burns holes, the 
attenuation will be less than aoL. Practically, the equa- 
tions will be applicable for aoL not too much greater 
than unity. 

Finally, we point out that dhole and Fhoio, as defined 
in Eq. ()14p . refer to the structure in the population, not 
to the depth and width which would be measured in an 
absorption experiment. From (|14p we always have the 
relation: 

1 - dhole = ■ (18) 

V i hole / 



2. Absorption and phase shift from a spectral hole 

Let us now employ Eqs. ([4]) and p?| to calculate the 
attenuation and phase shift of a laser field in the presence 
of a spectral hole. We take for the 5(A) function: 



The division by G (the fraction of atoms in |g) far off 
resonance) assures that 5(A) is correctly normalized to 
unity away from the spectral hole, i.e. ao is the absorp- 
tion length for a weak laser field in the absence of the 
spectral hole. 

Let us retrace our steps so far and underline the ap- 
proximations made. We have reached the two important 



5(A) = 



1 - 



i^ + A2 



i%i- + (A-Ao)2 



(19) 



i.e. we have a spectral hole burned at frequency Aq into 
an inhomogeneously broadened Lorentzian profile with 
width Finh centered at A = 0. Inserting this into Eq. ^ 
we find, with the help from the residue theorem. 
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ai{uj) = ao 



rhoio(rhoio+rh) 

4 



dhoU 



A2 

^0 



(Fhoic+rh 



A2 

^0 



(r...+r.)- + (Ao - u;y 
I 



ao 



ao 



rhoic(rhoic+rh) ? 

4 "'hole 



4 



(Ao-c.)2 

hole 



(r...+r.o- + (Ao - c.^ 



(20) 
(21) 



with Fhoic and di^o\c defined in Eq. (|14p . In the first terms 
we assumed that T^h 3> Th- The arrows indicate the 
limit when Finh — > oo, i.e. when we neglect the effect of 
the possibly very wide inhomogeneous background. Note, 
we had to include the inhomogeneous profile in this cal- 
culation since otherwise (?(A) — > 1 for A oo, making 
the integral in Eq. ^ divergent. Note, also that the fac- 
tors containing Finh are given by Eq. ^ and they may 
be replaced by the corresponding results for a Gaussian 
profile. 

If we compare Eqs. p7)) and ((20)) we see that in an 
absorption measurement with a weak field (not changing 
the populations further) the measured width and depth 
of the hole are related to Fhoic and dhoio by: 



of the complex Rabi frequency: 



p(moas) y. 

^ hole ^ ^ hole 



F. 



j(mcas) 
^holc 



FholoC^l 



hole 



Fhole + Fh 



(22) 



D. Calculation of error signals 

Now, let us turn to the calculation of real error signals 
used in the locking procedure. We will start simply by 
connecting the Rabi frequency, Q, used in the theoretical 
calculations, with real experimental parameters, and we 
will also state the well known results from Pound-Drever- 
Hall frequency modulation [l3, . 

To this end, we assume the laser field propagating 
into the locking system is perfectly monochromatic, with 
£:(in)(i) = Relfoe"'^"*} in the rotating frame. The elec- 
tric field is connected to the incoming intensity through 
the relation = |£olV2AtoC, which again is connected 
to the incoming power, P, and cross-sectional area A by 
j(in) ^ p(in)/A. Since we have defined the Rabi fre- 
quency as r2 = fj,£/h, where £ is the complex electric 
field inside the atomic sample, we can equivalently de- 
scribe our monochromatic incoming laser field in terms 



J 



niO,t) ^ noe-'"^"' withf^o = T\/^^- (23) 

h V nA 

Here the refractive index n takes into account the fact 
that the electric field is different inside and outside the 
atomic sample. 

Next, we phase modulate this incoming laser field at 
frequency lu^ such that: 



0(0, t) = J7ge-*['^«*+'"'^'"('^»*)l 



die )e 



(24) 



where m is the modulation index, which is assumed to be 
small in the second line, such that only the first sideband 
is significant. Jo and Ji are the Bessel functions of the 
first kind with m as the argument. 

It is now time to pass this modulated field through 
the atomic sample. We shall simply assume that for the 
carrier, upper sideband, and lower sideband, respectively, 
there are constant phase shifts of (pc, 4>+, and 0_ and 
constant transmission coefficients rjc, 77+ , and ?]_ of the 
amplitudes, such that the outgoing Rabi frequency can 
be written: 



-iAnt 



(25) 



These i] and (f> parameters may be calculated from 
Eqs. ([20]) and ([2T|) . When the above field hits a photo- 
detector, they will generate a photo-current proportional 
to the power p(°"*)(i) incident on the detector. This 
power is time-dependent and can be calculated from 
the above equation together with the absolute square of 
Eq. ^ to give: 



iv+ +f?-)sin (/«, 



Vc 

r 



+ 0- 



sin ^Wiiit -|- - 

cos oj^t + 



(26) 



It takes some trigonometric relations to reach this, and we have neglected the higher order contribution oscillat- 
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ing at 2u;n 



Harmonic analysis 



The first term in square brackets is a DC term de- 
scribing the average power reaching the detector. The 
second term in square brackets is that useful and impor- 
tant for laser stabilization and closely resembles the well 
known result from 's^. Here, this term is not expressed 
in the most useful form for calculations, but it clarifies 
some important qualitative properties of the error signal 
at different demodulation phases for our particular setup. 



First of all, the term 



is always rather small. 



For a broad inhomogeneous profile (Finh 3> Wm) in the 
absence of hole burning the phase changes little for a fre- 
quency variation of Wm, and the phase difference 0+ — (f)^ 
is small. In the presence of a spectral hole the phase may 
change dramatically with small frequency changes, but 
the two side holes at frequencies Aq ± Wm are almost 
identical and the phase difference remains small. In any 
case, we always have < aoLuJni/^inh <C 1 (un- 

less aoL is ridiculously large). As a result, we need only 
consider the two terms in the rightmost square brack- 
ets above to be pure sine and cosine terms, respectively. 
Note, some work has been performed on narrow inhomo- 
geneous profiles in the case when Finh ^ does not 
hold M. 



Secondly, the term < 



'^+^'^^ depends strongly on the 
situation. If the atomic contribution to the phase varies 
slowly versus frequency on a scale broader than cjni, this 
terms is small. For instance, for an inhomogeneously 
broadened profile with width Fjnh and no hole-burning 
effects, we must have 0c — '^^1^'^^ = 0([t^m/rinh]^), since 
4>c = "^^^"^^ to first order. In this case the sin(a;nii)-term 
in Eq. is second order in Wm/rinh- At the same time, 
we have in general that 77+ — ?]_ varies to first order in 



^m/rinh, while cos( 



't>++<P- 



) « 1. Hence the cos(cji„t) 



term is much larger than the sin(a;int) term, which can 
be neglected. 

The situation changes when hole burning is present. 
Then, in general, the difference between and 4>± can 
be quite large, and both cos((/)c — ^±±^) and sin((/)c — 
'^+~^'^- ) can be on the order of unity. At the same time, 
T7_l_ + r]- can be on the order of unity while 77+ — ry_ is of 
the order of Wm/rinh- Hence the sin(a;mi) term is on the 
order of unity and dominates the weaker cos{ujait) term. 

Experimentally, we may demodulate the signal given 
by Eq. (j26p by mixing it with a local oscillator in order 
to obtain a useful error signal for the laser stabilization 
system. When choosing the phase in this procedure there 
may be unknown phase shifts in amplifiers, etc., but the 
"cosine" or "sine" quadrature can be found by maximiz- 
ing or minimizing the error signal when hole burning is 
not present (for rare-earth-ion-doped crystals this can be 
done at an elevated temperature). 



Above we have calculated the detected power in a case 
when the input light is perfectly monochromatic. In re- 
ality the input light will vary in amplitude and frequency 
over time. The amplitude variations are relatively easy to 
measure and correct for; we briefly discuss how we do this 
experimentally in Sec. IIII Al Hence, in the following we 
concentrate on frequency errors. We may imagine that 
the transmission coefficients ijc, r]+, and ?7_ and phase 
shifts (jjc, (/)+, and 0_ in Eq. (|26p vary in time, which 
calls for further calculations. In general, a convenient 
method is to assume the incoming laser field to be of the 
form: 



(27) 



i.e. we have an almost single- frequency laser at Aq, but 
with an additional small harmonic disturbance of the 
phase with frequency uj and magnitude e, which we as- 
sume to be much less than unity. We note that the in- 
stantaneous laser frequency, Wg"^', is then given by: 



UJ, 



inst 



d<f>{t) 
dt 



Aq -|- ew cos{uut). 



(28) 



The above model is valid when the laser is running in 
the frequency stabilized mode with a narrow line-width. 
After phase modulating the beam at frequency Wm the 
light field becomes: 



f](0,t)«f]o(l + |e-'"*-|e'"') 



X (Jo + Jie" 



Jie''^-*)e 



(29) 



which is just Eq. (p4)) together with the harmonic error 
of Eq. (Hil) for e< 1. 

We now have nine frequency components which will 
propagate through our atomic medium, and we will treat 
these in a similar manner as in Eq. (|25p . To this end we 
make the following assumptions. Since e is small there 
is so little energy in the w-sidebands that these have no 
hole-burning effects. Consequently, we will have three 
spectral holes at frequencies Aq and Aq ± Wm- For these 
holes we have Fhoio ^ such that the holes can be 
treated independently when utilizing Eqs. (^0]) and (PT|) . 
We also neglect the effect of the inhomogeneous profile 
(Finh ^ Wm)- As a conscqueuce, the two spectral holes at 
frequencies Aq ± become identical and we introduce a 
common absorption coefficient, rys = 77+ = ry_ and phase 
shift, (f)s = (t>+ = (t>- where "s" refers to the sidebands at 
frequency cj,n. 

Introducing the hole widths rhoie,c, rhoio,s and hole 
depths dhole. c, c'hoio.s for the carrier and sidebands, re- 
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spectively, we may define: 



2. Low-frequency errors 



r].i{uj) = exp 



1 - 



Thole, i(rhoic.i+rh) J 

4 ^holci 

(r.„.e..+r.)^ ^^ — 



T Thole, i J , , 

a^L 2 "holc.iW 

2 (rho..,.+r,)f--- 



where i = c, s, and w now refers to the positive uj- 
sideband of the harmonic error, i.e. 770(1-^) is the ab- 
sorption of the frequency component Aq + oj. We note 
that rii{—ijj) = rjiiuj) and 4>i{—uj) = —(jjiiuj). Eqs. pO]) 
and pTjl have been derived directly from Eqs. pp]) 
and (I2ip with appropriate frequency substitutions. For 
instance, for the hole at Aq + we must substitute 
Ao Ao + Ldm in Eqs. and but, on the other 
hand, for the tj-sideband we must use Ao+ajm+w instead 
of Lo as argument in Eqs. (pO)) and (HI]). Consequently, 
Eqs. pO]) and ((3T|) depend only on cj. After a little alge- 
bra, the outgoing field can be written: 

Q.{L,t) = Jono[rjciO) - ie'nduj) sm{Lut + (t>c{uj))]e-"^"' 

+Jir!o[?7s(0) - zer/s(^) sm{u;t + 0,(tj))]e-*(^«+"»)* 

- Jir!o[?7s(0) - zer/s(^) sin(c^t -I- 0,(tj))]e"*(^«--'")*. 

(32) 

Assume that this field impinges on a photo-detector and 
let us calculate the power versus time, as we did in 
Eq. (|26p . The term oscillating at frequency ojai can be 
written: 

P!.THt) = 4F('")JoJisinK,t)x 

rv,e< : ■ ewe 



(33) 



where we have neglected terms of order O(e^). Firstly, 
we note that this is proportional to JqJi and oscillates as 
sin(a;mi)- In the example discussed around Eq. (^5]) we 
would have found similar behavior if we had assumed the 
side holes to be identical (with cj)^ — and ry+ = rj-). 
Secondly, in the curly brackets the real part of the factor 
ewe*"* is just the instantaneous frequency of the incoming 
laser (relative to Aq), see Eq. ([28]) . Hence, the fraction 
in the curly brackets acts as a transfer function map- 
ping this harmonic frequency excursion onto the mea- 
sured power (and eventually the output voltage) on the 
detector. This is similar to the way in which a complex 
impedance Z{u!) maps a complex current I{uj) onto a 
complex voltage V{u!) = Z{uj)I{uj) for individual Fourier 
components in electrical engineering. Hence, the trans- 
fer function T{u!) is directly applicable for purposes of 
feedback loop design for the laser stabilization system. 



Before we proceed to utilize Eq. ((33)) . we will make 
a small adjustment. We have assumed that the spectral 
holes burned at Aq and Aoiwm are static. This is a poor 
approximation if we consider harmonic errors which vary 
slower than the hole lifetime, T^g, so instead we will make 
a more general assumption that the actual hole center, 
Ao(t), is time-dependent, and is given by the history of 
the instantaneous frequency of Eq. ([28]) averaged back in 
time with a characteristic memory time, Tj-g: 

Ao(i) - 7^ /* wr*(t')e-(*-*'^/^-di' 

-'rg J-00 



= Ao + Re 



1 



1 + iivTig 



(34) 



We also assume that the hole shape is otherwise un- 
changed. This is only a good model for small frequency 
excursions. Now, instead of calculating the frequency 

Aq we will use: 



excursion from the hole center as ui™^*" 



Ao(i) = Rc 



ILUTr 



1 -|- zwTi 



- ewe'"* 



rg 



(35) 



This is the effective instantaneous frequency excursion 
from the hole center, which becomes small compared to 
Eq. (pS]) for low frequencies. To obtain the transmitted 
power, Pi™*^(t), the complex frequency excursion in the 
curly brackets of Eq. ((35]) can therefore replace the factor 
ewe*"* in Eq. ([33]) yielding: 



T(w) 



77e(w)77,(0)e**=(") - ?7c(0)77s(^)e**=(") 



iuj + 



P^™t)(i) = 4P('") JoJiRe{T(w) - ewe*"*} • sin(w„,i). 

(36) 

For high frequencies this change plays no role, and hence 
Eq. (j36p is a very useful model for the atomic response 
to harmonic errors in laser frequency on all timescales. 
However, the extra ad hoc term =^ in the denominator 
is far from giving the full picture of laser stability at low 
frequencies. This is discussed further in Sec. Ill El 



3. Evaluating the transfer function 



The transfer function in Eq. (|36p can be evaluated by 
using the expressions in Eqs. (|30| and (|3T]) . Below we 
will consider different regimes of this function. 

Firstly, we note that by inserting Eq. (fT8]) into (|30p the 
transmission coefficient at the hole centers can be written 
r]i{0) = exp(— p^^^' . ). Secondly, it is clear from the 

denominator of T(w) in Eq. (|36]) that the inverse hole 
lifetime, T^^, is one characteristic frequency, and from 



Eqs. dsn]) and ^ that the hole width, Fhoici, 
characteristic frequency. 



is another 
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FIG. 1: The function f{x) defined in Eq. (|38j- At x = 2.41 
it attains its maximum value of 0.172. Physically, f{x) is 
proportional to the slope of 4>i{u}) at w = in Eq. (jSTJ. 



Now, let us evaluate the transfer function in the regime 
uj <^ Thole, i- Here we may use the approximation 77,: (oj) w 
r]i{0) and: 



UJ 



1 h 

where the function / is defined by: 



Tho 



lc,i 



hole.s 



!(a; + 1)' 



(37) 



(38) 



This function /(x) is shown in Fig. [T] If cj ^ rhoio,?: we 
will always have (j)i{uj) <C chqL, which means practically 
that the exponentials in Eq. p6p can be approximated as 
1 + i(t>i{Lo). We then find the transfer function: 



^(- + -)[/(^c)-/(:r,)]. 



In the central regime, T^„^ ^ w <gc F 



Fh 1 + iwTrg 

(39) 

hoie,i, the term 

, , . „ is unity and the transfer function is real. In the 

low- frequency limit, oj <^ 2^rg^' the transfer function be- 
comes imaginary and proportional to lo. 

In the high-frequency limit, lo 3> Fhoic, we have 



■qi{uj) = exp(- 



QqL 



), 01 (^^) 



4tij 



UJ 3> — . Then we obtain: 



T{lu) k. e 



<C 1, and 



(40) 



The transfer function T{uj) has been plotted in Fig.[2]for 
a choice of reasonable experimental parameters (see the 
figure caption). It is clear that there are three distinct 
regimes, as calculated above. Assuming that the terms 
in the square brackets in Eqs. ((39)) and ([40]) are posi- 
tive, the transfer function is a negative real number times 
iu), 1, and J- for the low-, medium-, and high-frequency 
regimes, respectively. This behavior is clearly seen in the 



10" 



10" 




T3 
0! 




10' 



102 lO'* 

Frequency (Hz) 



10*' 



FIG. 2: The magnitude and phase of the transfer func- 
tion T{ijj) calculated from Eq. (|36p (heavy lines). The light 
lines show the asymptotic cases discussed around Eqs. p9|) 
and (|40p . Parameters used are Ti = 150 \x.s, T2 — 18 
\LS, Trg = Tgr — 4 ms, 6cr ~ 0.5, flo — 2n -1 kHz, and 
m = 0.40, giving Fh = 27r-17.5 kHz, Fhoicc ~ 27r-21 kHz 
and Fhoica ~ 27r • 18 kHz. These parameters are close to use- 
ful working values, as we shall see in Sec. IIVI 



magnitude of T{lu) shown on the upper plot in Fig. [21 
The fact that the transfer function is real at medium 
frequencies means that the error signal (x Re{r(ti;)e*"*} 
will oscillate in phase with the actual frequency error 
oc Reje*'^*}. For high frequencies, the extra i factor 
makes the error sig nal oscillate as Re{e'['^*-''/2l}, i.e. the 
response is 90° delayed. This is shown as the phase reach- 
ing —90° in the lower plot in Fig. [21 For low frequencies 
the situation is the opposite; the phase is advanced by 
90°. This behavior of the gain and phase has been previ- 
ously reported in experiments and numerical simulations 



In our calculations we always assume that the power in 
the carrier beam is higher than in either of the sidebands, 
leading to Fhoic,c > Thoics- Then, according to the defi- 
nition in Eq. (|38p , Xc will be larger than Xs , and the term 
in the square brackets in Eq. will be positive, as we 
assumed above. For Eq. p9p . however, we can have a 
situation where f{xc) < f{xs) if, e.g. 2.41 < Xs < Xc, ac- 
cording to Fig. [Tl In this case there is a 270° phase shift 
between the medium- and high-frequency regimes which 
in practice means that the sign of the error signal can- 
not be chosen correctly for all frequency components in a 
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closed feedback loop. Physically, the sign change occurs 
when the slope of (j)c{^) around lj — Q equals the slope of 
(/)s(a;) in Eq. ([31]) and we must assure that this is never 
the case. Note also that prior to and in the initialization 
of the laser locking feedback loop the spectral holes are 
broad and shallow (since the laser is broadband). The 
carrier hole will be deeper than the sideband holes while 
the widths are roughly the same, limited by the broad 
laser line-width. This in turn assures that the slope of 
4>{iS) for the carrier is larger than for the sidebands. We 
must choose the right parameters such that the sign will 
remain correct when the feedback loop is closed and the 
laser line width narrows. 

We should also comment on our theoretical approxi- 
mation that e < 1 in Eq. According to Eq. ([28]) 
this means that the instantaneous frequency excursion 
at frequency lo must be much smaller than euj. At first 
sight, this seems to limit the usefulness of the calculations 
at low frequencies. However, for lo < Fhoic the linearity 
of ([?T|) versus cu enables us to relax this condition on e. 
In practice, the theory is valid when euj < Fhoic- 



4- Parameter choice in general 

The theoretical observations above enable us to dis- 
cuss the optimum parameters in general. Further estima- 
tions associated with our particular setup will be given 
in Sec. Ell 

We start by focusing on how to obtain a large error 
signal. The main reason for this is to obtain a high 
signal-to-noise ratio in the error signal. Our first ob- 
servation is the fact that the detected power given in 
Eq. is proportional to P^'"). It is no surprise that 
more light gives a higher signal at the detector, but it is 
wrong to just naively increase the incoming light power 
P*^'"^ and expect a better performance. Doing so will in- 
crease So in Eq- (|15p and in turn the hole widths Fhoic, i in 
Eq. (|14p . However, increasing the incoming power and 
at the same time increasing the beam area. A, leading 
to an unchanged intensity, will always help. Hence, it is 
a good idea to use an atomic sample with a large area 
orthogonal to the direction of beam propagation. 

Next we observe in Eq. ([55)) that the front factor JqJi 
attains its maximum value of 0.339 when the modulation 
index is m = 1.08. This value is often used in laser sta- 
bilization setups utilizing optical cavities [s^ . However, 
as opposed to the resonance lines in a cavity, the shape 
of spectral holes depends on the optical power. If the 
carrier and sideband powers were equal the holes would 
be identical, i.e. rjc = rjs and (j)c = (t>s, leading to a zero 
error signal according to Eq. ([55)1 . For this reason the 
optimum modulation index is lower than 1.08, leading to 
a more asymmetric power distribution between the car- 
rier and the sidebands. We have searched our parameter 
space with numerical methods while optimizing the sig- 
nal in Eq. ([55)1 . The result is that m = 0.56 is a good 
choice. However, there are further complications regard- 



ing the laser stability which suggest that m should be 
even lower. This will be discussed in Sec. Ill El 

Regarding the hole widths Fhoie,c and Fhoio.s it is clear 
from Eq. pQ)) and Fig. [1] that Xc = Fhoie.c/Fh should not 
be much greater than 2.41, since a higher value simply 
makes the /-function decrease again. Also, we would 
like to make Xg — Fhoic, s/Fh small in order to decrease 
f{xs). We can do this by lowering the modulation index 
TO. Making to too small will also lower the factor JqJi 
and this is why we found to = 0.56 to be the optimum 
choice seen solely from the point of view of optimizing the 
error signal. However, the magnitude of the error signal 
is not everything. The narrower the hole widths Fhoic, i, 
the longer the duration of the atomic phase memory and 
hence potentially better phase stability of the laser can 
be obtained. We should also note that a given width, 
Fhoic, i, can be obtained in different ways according to 
Eq. ([13]). One could choose a high intensity (high sq) 
and a short hole lifetime Trg (low R according to ([H])) 
if adjustable. On the other hand, a low intensity and a 
long hole lifetime could give the same result. In general, 
the latter will give the better long-term stability of the 
spectral hole. 

To estimate the optimum optical density aoL let us 
assume that Xc ~ 2 and Xg ~ I. This is not far from 
optimum given the discussion above. Inserting this into 
either Eq. ([5^ or ([^D)) leads to the ballpark estimate 
aoL w 1.15, corresponding to a background intensity 
transmission of e~"°^ « 32%. Note, this is on the edge 
of our approximation that agL should not be too large 
for quantitatively correct results. 



E. Laser drift 



In the previous sections we have calculated the error 
signals for laser locking based on the linearized model 
with the time-independent distribution function 5(A). 
For slowly varying errors on timescales slower than the 
hole lifetime, Tj-g, we presented in Sec. Ill D 21 an ad hoc 
model to account for the loss of gain. However, this does 
not really illustrate the real challenges in long-term sta- 
bility of the laser frequency. In the present section we will 
show that under certain conditions there is a solution to 
the equations where the laser is locked, but the frequency 
is drifting linearly with time. Below in Sec. Ill ETl we will 
re-calculate the shape function g{A) and then the ab- 
sorption lengths apt and ai in the presence of laser drift. 
This requires that we reconsider Eqs. ([Sl llip in more de- 
tail. Based on these calculations, we will in Sec. HIE 21 
derive corrections to the error signal from laser drift. 
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1. The drift model 

We consider a situation where the incoming laser field 
is given by 



The change in sign arises from the change of variables, 
and the saturation parameter now depends on A as: 



Ao +I3t, (41) 



i.e. the instantaneous frequency Winst varies linearly with 
time, where (3 is the drift rate in rad/s^. We assume that 
the drift is much less than Fh during an optical coherence 
time T2 (i.e. /3 <C F^). Then the populations pc and pg 
can be safely taken as constants on the timescale of T2, 
and the coherences u and v will follow the incoming field 
adiabatically. We find: 



u — IV 



-^»oe-'(^°+4^)*(pe-pg) 
^ + i[A - Ao - I3t] 



^ h 



(42) 



Inserting this into Eqs. (Pl fTTj) gives the population equa- 
tions: 



siA) 



p2 



^ + (A-Ao)= 



(50) 



The above equations describe the populations when the 
laser has instantaneous frequency Aq while the drift rate 
is p. The reason for multiplying by ^ becomes apparent 

when we now define the series expansion pc = p^'^ + 
^Po^'' +i^^/Oc^' + . . ., and similarly for pg and p^. Inserting 
these into Eqs. (|47ll49p leads to iterative equations for the 
different expansion coefficients: 



dpg 1 fs{t). 1 . A 1 1 



dpr bcr 1 1 



(44) 
(45) 



'dpi^^/dA 
dpi'^/dA 




Thole 



A(A) 



■u+iy 

Pe 

.0+1) 



Pr 



(J + I) 



(51) 



where s{t) is the time-dependent saturation parameter 



s{t) 



^ + (A - Ao - pty 



(46) 



where the dimensionless matrix A is given by: 



and So is still given by Eq. (|T5|) . Of course, the above 
equations can be solved by numerical methods. However, 
our aim is to derive an intuitive condition for the pres- 
ence of a linear laser drift. This is best done analytically. 
First, we note that the population equations (|43ll45p de- 
pend only on t and A in the combination of A — (3t. We 
will change variables to A. We will also make a series ex- 



pansion in the dimensionless parameter ^ 



This 



parameter is a measure of how far the laser drifts during 
a hole lifetime, T^g, compared to the width of the hole, 
Fhoic- We multiply the population equations by ^ and 
find: 



^dpe 
^dA 

^ 9pg 
^dA 



dpr_ 

^dA 



Lrg 



T'lFholc 

T.g 



s{A) 



TlFholo 

T,g Pg 

T^gr Thole 



2 

,s(A) 



[Pc -Pg]+. 

[Pc - Pg] + bcgPc 



TiF 



hole 



Pe 



Thole 
Trg Pg 



Thole 



(48) 
(49) 



'Trg (s 
Ti V2 



1 



r, 2 " 
_2k 1 

1 1 



(52) 



We have omitted A in s and A for brevity. Note, that 
the first and second rows in this matrix refiect Eqs. ([47|l 
and (Uni), respectively, for the individual terms in the 
series expansion. A solution fulfilling these two equations 
also satisfies Eq. (pS)) automatically. The third row in A 



(3) 



(47) 

ensures that p, 
the condition p, 



(0) 



■^^+pi>^ 



P^' 



Oforj > 1. Together with 
1 the total population 



is always conserved independent of ^. 



In order to find p^^ , pg^'' , and pr^'' , we need to calculate 
the inverse of A and we must know the zeroth order 



.(1) 



steady-state solutions po°^ , Pg°^ , and p^"^ . 
defined as in Eq. we derive 
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JO) 



G 



Ti 



\{1 + R) + 1T,, 



s 

'-^ • 2 

Ki + i?) + i' 



1 + ^ 



-1 - 



Ti 



G-\ 



Ti 2 
Ti V2 ^ 



1) 



^(s + 1) 

1] 



Ki + i?) + i' 



(0) 



Ti 



Ti 2 

Ti ^2 
[^(§ + 1) + ^"'^ 



1) 



Ti 2 J 

T„6er . 



Ki + i?) + i 



(53) 



(54) 



Here Eq. (f54|) is simply the steady-state solution of 
Eqs. pfTT|l or solutions of Eqs. (|47ll49p when ^ = 0. We 
also note that the zeroth order populations can be found 
simply by taking 



Pe 




0' 


Pg 


= A-1 





(0) 




1 









Taking the derivative of pe and pr with respect to A 



in Eq. (IM)) and multiplying Eq. ([5T|) from the left by 

FhoicA"^ we obtain expressions for p^c \ Pg^\ and pi^\ 
We remember that the distribution function is given by 

.9(A) = i(pg - Pe) = ^(pi°^ - pi°^) + ^(pi^^ - pi'^) + 

0(^^). The term linear in ^ is our correction, Adrift (A), 
to the shape function g{A) due to the drift. Carrying 
out the calculation we find: 



gdrift(A) = 



(3T^ sor2(A-Ao)[^ + (A-Ao)2] 



2(1 + 4^^)2 



ai(Ao) 



- holi 
4 



■ + (A-Ao)2 



1 3 



3 Tig 
2 T 



1 T2 
^ ig 

2T1 



T 

^rg 

2ri 



C?holo X 



Ti 


1 1 -^rg _|_ i rg 

^ 2 Tgr 2 


+ 4 


6er(l + ^) 


^rg" 

''<=gTg,_ 











(56) 



(57) 



where Fhoic and dhoio are defined in Eq. ()14|) . The second 
line is derived with help from the residue theorem by in- 
serting Adrift (A) into Eq. Q and setting w = Aq, i.e. we 
have calculated the imaginary absorption length ai ex- 
perienced by the laser beam at frequency Ao drifting at 
rate /3. The first-order drift contribution to Q!r (Ao) is 
zero, which is easily seen by symmetry. 



2. Error signal from drift 

Expression ([57]) may look a little complicated, but all 
we really need to understand is the factor to the left of 
the x-sign. This factor says, that ai should be calculated 
by taking the background ao times "how far we climbed 
up the hole" (this is ^ = pp^) times the relative depth, 
echoic, of the hole. The rightmost fraction in Eq. ([F7|) 
is merely a constant independent of laser power. Hence, 
this constant is the same for the center and side holes. In 
the, not so uncommon, case when focr^ig ^ Ti (leading 
to i? ^ 1) the factor is approximately equal to [2(1 + 



Y^)] ^- For brevity we make this approximation below 
remembering we can always return to the more accurate 
expression if necessary. 

In order to calculate the error signal from a drifting 
laser, we now add a phase modulation at frequency 
with modulation index m to the incoming field of Eq. (j41[) 
in a completely analogous way to Eq. ([24]) for a static 
incoming field. This again leads to a carrier and two 
sideband fields and we may employ Eq. to calculate 
the error signal. When the drift is slow (^ <C 1) the 
absorption coefficients r]c, 77+ and r]_ will be given by 
m — exp(— ^77^ p ), and the phase shift will be given 

by times Eq. ([57]) . We then obtain: 
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1 + ^ 



^holc,c 



- holc.i 



^holc,s 
hole.s 



[■holes J 



(58) 



where the approximation sin(0c — 0s) ~ 0c ~ 0s is em- 
ployed. Comparing this expression to the low-frequency 
version in Eq. p9p we find most importantly that the dif- 
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0.4 0.6 
modulation index m 



FIG. 3: Thresholds for stable laser operation. If d'^^^l'' is 
below the solid line the condition in Eq. (|59p is met and the 
zero-drift solution is stable. In order for the low-frequency 
part of the laser locking to have the correct sign {f{xc) — 
f{xs) > in Eq. (I39|l ') we require the less stringent condition 



i(meas) 
•'hole.c 

the laser can be locked in a linearly drifting mode. 



that rfj^oij, J is below the dotted line. Between the two lines 



ference in /-functions has been replaced by the difference 
in the ratios dhoic.i/rhoic.i- As discussed previously, there 
is a risk of obtaining the wrong sign for the error signal. 
The difference in square brackets in Eq. (|58p must be pos- 
itive for zero drift with /3 = to be a stable solution. If 
this is not the case, there will be a bi-stable solution with 
positive or negative non-zero values of /3. To calculate the 
magnitude |/3| of this drift rate requires a complete cal- 
culation of the non-linear problem of Eqs. (|43ll46p which 
is outside the scope of this work. 

We wish to operate the laser stabilization system with- 
out drift, and to this end we derive a stability criterion 
based on convenient experimental parameters. The drift 
problems occur when the hole depths are too large, as 
we shall see below. Given the modulation index, m, we 
will determine a threshold value that the measured center 
hole depth, d^'^^^l^ , should not exceed. The hole width 
and depth can be measured by scanning a laser across 
the hole and monitoring the absorption, as we will show 
in Sec. IIVI Also, for practical values of the modulation 
index, the carrier has the major fraction of the optical 
power, and the total light transmission of the locking 
beam itself can be utilized to obtain an approximate es- 



timate of d 



(meas) 
holc,( 



We first observe that r^^j^, ^ 



^ h 



where we define y = Jf/Jg as the power ratio be- 
tween one of the first sidebands and the carrier beam. 
To see this, use Eq. remembering that sq has to 
be multiplied by Jq and Jf to calculate the correct 
saturation parameter for the carrier and sidebands, re- 
spectively. In addition, with Eq. and the thresh- 
old condition dhoio,c/rhoie,c > rfhoic.s/rhoie,s, we derive 
(rhoic,s/rhoic,c)^ > y^^^- Next, we replace Fg^j^ .. by 
yi^Y^oic c ~ -I- r^, and we are left with a condition de- 
pending on only the two parameters y and Fhoic,c/rh. A 



combination of Eqs. ((18)) and (1221) will show that the mea- 
sured hole depth can be written 6?^™!^'*^^ = 1 — Eh/Fiioicc 
and we finally reach our criterion: 



,(meas) , 
"hole.c ^ 




y = 4- (59) 

•^0 



Given the modulation index, m, we can calculate y nu- 
merically and plot the threshold value as a function of 
TO. This has been done in Fig. [3] (solid line). 

We may also derive a criterion for ensuring the cor- 
rect positive sign of the factor [f{xc) — f{xs)] in Eq. ([39]) 
assuring that the low frequency part of the error signal 
is correct. There is no simple analytical expression for 
this, but a numerical result is shown by the dotted line 
in Fig. [3l Note, that it is easier to fulfill the criterion 
for correct low-frequency behavior than the criterion for 
no linear drift. This is an important observation which 
shows that all the calculations regarding the drift model 
are worthwhile and necessary to obtain a complete un- 
derstanding of laser stability. It is indeed possible that 
is in between the dotted and solid lines in Fig. [31 



,(mcas) 
''hole,. 



in which case the laser stabilization system is apparently 
locked but still the laser is drifting linearly. 

We conclude this section by pointing out that the drift 
calculations can be performed in a similar manner for a 
simple two- level system in absence of a reservoir state |r) , 
but the results can readily be guessed by setting Tig = 0, 



oo, be 



1, and &p 



0. Then the three-level case 



will reduce to the two-level case and Eq. (|57p will reduce 
to ai(Ao) = -^j^dhoic 



F. General remarks on the calculations 

Up until now we have considered a two-level sys- 
tem with a single reservoir state to model the trapping 
of atoms in the hole-burning process. This is a sim- 
ple system which allows for not too complicated ana- 
lytical solutions, thereby maintaining the physical un- 
derstanding. This simple system is actually found in 
Tm'^+:Y3Al50i2 [H, [2§|, and the even simpler pure two- 
level system is found in Er^+iYaSiOs and Er:KTP 

However, more complicated cases exist. For our exper- 
iments with Pr'^+:Y2Si05 there are three ground state 
levels and three excited state levels, see Fig.[ll[e), and all 
nine possible transitions exist, due to the inhomogeneous 
broadening, and in principle play a role. If one wishes 
to extend the theory to cover this we must in Eq. (Q 
replace aog{^) by ^jaojgj{A), where j identifies the 
individual transition, gj{A) is the shape function calcu- 
lated with the appropriate saturation parameter given 
the transition strength, and aoj is the absorption length 
for the individual atomic species. The total absorption 
length must fulfill ao = J2j o^Oj- 

To calculate 5j(A) for the case with three distinct 
ground states we consider the physical systems shown in 
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a) ±1 b) _l£). 




FIG. 4: Different level schemes used in this paper. We de- 
fine timescales for relaxation between ground state levels and 
branching ratios from the excited state. The excited state 
lifetime is always denoted Ti. (a) The most naive scenario 
with two levels, considered in Sec. Ill Bl and in the first row 
of Tab.|T] (b) Our basic model for all calculations, described 
in Sees. Ill CI and III El Rows two to four in Tab. |I] refer to 
this case. Note, we may have different relaxation timescales 
Trg 7^ Tgr. (c) and (d) Different schemes with three ground 
states coupled as shown with RF-magnetic fields. Hence the 
timescale is the same in two opposite directions. These cases 
are reflected by rows five and six in Tab. IT] respectively, (e) 
The real Pr^+:Y2Si05 level scheme. 



Fig.mjCjd). The active optical transition is still |g) — + |e), 
but two reservoir states are present. In Fig.[5Jc) an asym- 
metric case is shown where the reservoir states are labeled 
|c) and |f) for "close" and "far", respectively, describing 
their position in the RF pumping scheme relative to the 
state \g). In Fig.l^d) the symmetric case is shown where 
the reservoir states are labeled |rl) and |r2). It is a sim- 
ple matter to generalize the methods of Sec. Ill Ci to three 
ground state levels. Doing so, we calculate steady-state 
populations in order to derive the distribution function 
(7j(A). This function is found to have exactly the same 
form as Eq. apart from new values of R and G, 

which are given in Tab. U (fifth and sixth rows). As shown 
in Fig. |3tc,d) we assume the same timescale in both di- 
rections between two ground states, e.g. Td — Tfc, etc. 
This is valid in our experimental case with RF-magnetic 
fields coupling the adjacent levels. 

Hence, for a single atomic species, adding more ground 
reservoir states only changes the spectral holes quantita- 
tively, but qualitatively we still have a Lorentzian-shaped 
hole fulfilling Eq. (|14p. as for the simple case of two- levels 
plus a single reservoir state. However, for multiple atomic 
species the contribution of different transitions with dif- 
ferent strengths may lead us to sum up quite different 
f/j(A) functions with a non-Lorentzian shape as the re- 
sult. Instead of performing a complete quantitative ex- 
amination of this, we tried in the experiments to keep 



Case 


G 


R 


_L W yj ic V ci 


C — 1 


R — 1 


Three-level 




— 4t- 


Three-level, 
natural decay 


G = 1 


- 1 -I- ''"J^e 
J 1 


Three-level, 
RF eraser 


G 2 


R- 1 


Four-level, 
RF eraser (1) 


G = 1 




Four-level, 
RF eraser (2) 


G = 1 


R- T.l''o,l+T,.i>„,] 



TABLE I: The value of G and R for the different setups shown 
in Fig. m The first row gives the relations for the two-level 
atom (Fig. |4ja)) and the second row describes the two- level 
plus reservoir state system considered in Sec. Ill C] (Fig. l4tb) ) . 
The third and fourth rows are special cases of the second row. 
In the third row we assume Tgr = oo, which describes a one- 
way natural decay from states |r) to |g). In the fourth row 
we assume Trg = Tgr which describes the case when an RF 
magnetic field couples the otherwise uncoupled states |r) and 
|g). The fifth and sixth rows correspond to the cases shown 
in Fig. |4jc) and Fig. |4jd) , respectively, where there are two 
reservoir states. These four-level cases are presented since 
they resemble our experimental case using Fr^+iYaSiOgas the 
atomic medium. 



this difference small, in order to mimic the three-level 
system and demonstrate the qualitative features of the 
theoretical calculations. The results in Tab. U will help 
us do this. 

We also wish to remind the reader that our theory gen- 
erally assumes perfect lasers or perfect lasers with har- 
monic errors. In real life this is not the case, but our 
approximations are still quite good if the stabilization 
system maintains a narrow line-width in the laser. If the 
laser line width, for example, is 1 kHz and the hole width 
is 20 kHz, there will be some kind of folding effect of the 
order of 5%. Also, we have assumed that aoL is not too 
large. If we, for example, set a^L « 1 and assume a mea- 
sured hole depth of around 50%, the transmission of the 
carrier beam is e^^/^ « 60%. This means that the satu- 
ration parameter varies by 40% over the sample, and we 
can approximately take this into account by lowering the 
saturation parameter to 80% of the value calculated from 
Eqs. ([15]) and (|23|) . In this manner (for the two levels plus 
a single reservoir state) we should be able to keep the the- 
ory quantitatively correct within around 10%, while all 
the qualitative features should hold true. 



III. SYSTEM DESIGN 

Based on our theoretical calculations we are now able 
to design the electronic feedback system giving the best 
performance. Below in Sec. IIII Al we describe our setup 
for the feedback system, and in Sec. IIII Bl we describe 
how to estimate the limitations of the performance. 
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FIG. 5: (a) Optical and electronic design schematics. Ab- 
breviations: N, notch filter; LP, low-pass diplexer; PS, phase 
shifter; LO, local oscillator; EOM, electro-optical modulator, 
(b) Electronic schematics of the RF-eraser circuit. See the 
text for component values and details. 



A. Experimental feedback system setup 

The main building blocks of the feedback system are 
shown in Fig. ^a,) . Starting from the upper right corner 
we show the laser cavity of our Coherent CR699-21 dye 
laser. The most important elements controlling the laser 
frequency in the commercial version of the laser are a 
piezo-mounted mirror and a Brewster plate. These were 
retained in our setup. In addition we placed a Linos 
PM25 electro-optical modulator (EOM 1) inside the cav- 
ity. We fed the two electrodes by two separate amplifier 
circuits; ICl which is fast and IC2 which can deliver high 
voltage (±200 V). 

From the laser output the laser beam was directed 
through a New Focus model 4002 modulator (EOM 2) 
applying uj^ — 27r • 50 MHz modulation from a local 
oscillator (Wenzel Associates, Inc.). The modulated 
beam was then expanded to cover the entire area of a 
Pr'^+:Y2Si05 crystal with diameter of 19 mm, thickness 
of 5 mm, and doping concentration 0.005%, purchased 
from Scientific Materials Corp. This crystal was kept in 
a cryostat (Oxford Optistat CF-V) operated at 3.0 K. 

The electronic circuit design for the RF-eraser coil sys- 
tem, which enables us to control the hole lifetime, is 
shown in detail in Fig. [51(b). Surrounding the crystal are 
two home built sets of coils denoted L2 and L3, each con- 
sisting of two three-turn coils with a diameter of 23 mm 
and 7 mm separation. These coils function practically 
as a transformer, coupling the 10.19 MHz and the 17.31 
MHz channel, since they are placed on top of each other. 



In order to effectively reinforce inductance, making res- 
onant circuits out of the two RF channels a band-stop 
filter consisting of Li = 0.76 |i.H and C5 = 100 pF is 
inserted and adjusted to block the 17.31 MHz signal. By 
adjusting Ci and C2 in combination we can now tune the 
resonance frequency for the 10.19 MHz channel and as- 
sure a 50 n impedance match. Likewise, C3 and C4 are 
adjusted to assure the correct resonance frequency and 
impedance of the 17.31 MHz channel (the two channels 
are not completely independent and adjustments must be 
made iteratively). The effective values of Ci - C4 depend 
to a large extent on stray capacitances in the approx. 50 
cm leads from the outside adjustable capacitors to the 
coils inside the crystal. We used 35 pF adjustable capac- 
itors for Ci - C4 in addition to 22 pF, 22 pF, and 47 pF 
capacitors coupled in parallel to Ci - C3, respectively. 

The RF fields are generated as sawtooth sweeps, the 
10.19 MHz signal is 100 kHz wide, and the 17.31 MHz sig- 
nal is 200 kHz wide. The sweep time is 0.82 ms, which is 
comparable to the hyperfine level coherence time of 0.50 
ms and hence the pumping becomes effectively incoher- 
ent. This procedure assures smooth re-population over 
time of the hyperfine levels since atoms with different fre- 
quencies on the inhomogeneous hyperfine transition are 
affected at different times. 

After the light has passed the crystal it is measured 
by a Thorlabs detector (PDAlOBS-AC, bandwidth 150 
MHz) and the output is sent through a notch filter to 
block possible higher order modulation at 100 MHz, 
which is known to cause systematic offsets in cavity lock- 
ing setups ^ . The notch filter also blocks possible elec- 
tronic pick-up at 10.19 MHz and 17.31 MHz. The signal 
is then mixed with the phase-shifted local oscillator on 
a Mini-circuits LPD-2 mixer and after low-pass filtering 
(20 MHz cut-off), the error signal is produced at DC 
frequency. Both the notch and low-pass filters are de- 
signed as diplexers with 50 Q impedance matching at all 
frequencies to avoid reflections (this is particularly im- 
portant for the low-pass filter to avoid back-reflections 
into the mixer [42|). 

Based on the error signal we actuate the laser fre- 
quency using EOM 1 which is driven by the pure ana- 
log electronics shown in Fig. [Ha) around ICl and IC2. 
The complex electronic gain from the error signal to the 
voltage across the electrodes of EOM 1 is given by: 



5(^) 



R2 
Ri 



iuj{R3 + Ri)C + 1 



(60) 



We see that there is a characteristic cutoff frequency, 
fc — ^ = 2Ti-(fl34-_R4)c ' 'which separates this gain into 
a low-frequency part proportional to J- and a high- 
frequency part where the gain is constant and real. The 
resistors Ri to i?4 are adjustable but designed to work 
in a range such that the critical frequency, ujc, can be 
set equal to Fhoie- In this manner the electronic gain of 
Eq. (I60p together with the medium- and high-frequency 
part of the atomic response shown in Fig.[5]Ja) add up to a 
total response proportional to J-. The actual component 



16 



values used were Ri = 50 ft, R2 = 50 kil, R3 = R4 = 1 
kfi, and C = 6 uF. The resistors R2 - R4 are adjustable. 
The amplifier circuits are based on ICl, which is Ana- 
log Devices amplifier, model AD8021 and IC2 which is 
an Apex model PA85 amplifier. Circuits for maintaining 
low offset levels, for preventing integrator windup [43|, 
and for signal monitoring, are not shown in Fig.[5]Ja). 

The output of IC2 is also sent to a low-frequency unit 
which consists of a PC processing the input digitally and 
in turn sending a signal to the commercial part of the 
laser. Specifically, the input from IC2 is sent to a PI 
regulator (i.e. an amplifier with gain g{uj) oc |^ -I- C2 oc 
i/T+i(jj y chose T to be equal to a typical value 

for the hole lifetime, Tj-g. Combining this particular re- 
sponse of the digital processing with the atomic response, 
corresponds mathematically to replacing the denomina- 
tor ioj + J- in the transfer function in Eq. (1551) by ioj. 
In other words, the digital processing together with the 
low- and medium-frequency response shown in Fig. [^^a) 
amounts to a constant gain response. Together with the 
low-frequency response, proportional to of the analog 
electronic part, the total response of the entire analog and 
digital system at all frequencies becomesproportional to 
i . This is known from feedback theory [43] to assure sta- 
ble operation. The digitally processed signal is added as 
a current to the reference photo-diode behind the stabi- 
lization reference cavity. In this manner the commercial 
laser control system will detect a frequency error and ac- 
tuate the piezo-mounted mirror and the Brewster plate. 
Finally, to assure that the working point of the refer- 
ence cavity remains within an appropriate range, an in- 
tegrated version of the digitally processed signal is sent to 
the external scan input of the commercial laser control, 
which in turn changes the reference cavity length and 
maintains the current added to the photo-diode around 
zero on average. 

Apart from a stable frequency it is also important to 
stabilize the laser output power. To this end we split 
off a small fraction of the laser beam to monitor the 
laser power. This in turn is used in a feedback loop 
to control the amount of pump light sent into the laser 
gain medium. Our dye laser is pumped by a Coherent 
Verdi- V6 laser, and a Brimrose AOM (model FQM-80- 
1-.532/WQ) is placed in the pump beam to adjust the 
power. Details of the feedback electronics are available 
in 

B. Factors limiting the laser performance 

Two important design parameters in the setup are the 
group delay of signals in the feedback loop and the signal- 
to-noise ratio (SNR). For this reason we chose high-speed, 
low- noise amplifiers, and the group delay of the filter 
diplexers was carefully designed so as to be constant over 
relevant frequency ranges in order to avoid distortion. 
The total feedback loop delay time is below 100 ns. For 
technical details see [Sfl, IMj. Note also, that short de- 



lay times are only possible if the modulation frequency 
is high. 

As discussed in Sec. IIII A[ the total response of the 
atoms and the electronic feedback system is a gain func- 
tion proportional i . This is purely imaginary, and corre- 
sponds to a constant phase shift of —90° for all frequen- 
cies. In reality, there are corrections to this since am- 
plifiers, mixers, and filters have finite bandwidths. The 
total group delay, r, of the entire system can be mod- 
eled as a gain function, (7dciay(w) = e"*"^"^. If the total 
feedback gain is too strong we expect oscillations at a 
frequency where the feedback phase is —180°. This will 
occur around a frequency for which e^''^'^ — ~i since we 
then add —90° to the high-frequency phase shift which is 
already —90°. For r = 100 ns we estimate a critical fre- 
quency for oscillations /osc = ^^5^ = ^ = 2.5 MHz which 
agrees well with our experimental observations. We dis- 
cuss the consequences of this delay on laser stability in 
Sec. UVCl 

Let us describe how to estimate the impact of de- 
tector noise. For our Thorlabs PDAlOBS-AC detector 
the electronics noise specification is —63 dBm in a 1 
MHz bandwidth around the 50 MHz center frequency, 
which can be translated into a noise in the output volt- 
age of 0.16 -^S=. The shot noise current measured in 

V Hz 

a bandwidth B (in Hertz) is given by isn = \/2eBid 
(43 |. where id is the detected current. Knowing the 
power-to-current conversion, Cj/p — 0.37 A/W, and 
the trans-impedance of the detector, Cu/j = 1.65-10'' 
V/A, we may calculate the output voltage noise given 
the detected power, P. Adding this to the electronics 
noise we obtain U^°l'"' = 0.18 ^ v/P[mW] + 0.77. Af- 
ter the photo-detector we mix the output signal with 
the 50 MHz local oscillator and low-pass filter the re- 
sult with a bandwidth B = 20 MHz. We can thus 
calculate the equivalent noise power given the detected 
DC power, P as PS°r = US°r^/iCi/pCu/i) = 
1.3 • 10"''mWv/P[mW] + 0.77. This should be compared 
with the amplitude, Pf^f""^ = 4PJo JiRe {T(w)etje*'^*}, 
of the 50 MHz modulations in the power from Eq. ([36|) . 

As we shall see, in the experiments described in Sec. lIVI 
the best laser performance was observed using a detected 
power, P, of 0.23 mW which corresponds to the equiv- 
alent noise power = 1.3 • 10"'' mW = 5.7 • lO'^P. 
The modulation index, m, is 0.20, giving JqJi — 0.1. 
In the high frequency regime the transfer function T(uj) 
can be estimated from Eq. (|40|) . Assuming aoL = 0.66 
(corresponding to 52% transmission) and p^^^^" — « 0.5, 

■p^i- « 1, we find T(uj) w _M23. The signal can then 

Thole, s ^ ^ lU! ^ 

be estimated as P^^t"'*' ~ 6.4 • 10~'*P for a harmonic error 
with phase variations of one degree, i.e. e — 1 •27r/360. 
This estimate is close to the noise limit, which means that 
with a detection bandwidth of 20 MHz we can essentially 
detect a 1 degree phase error with a SNR of unity after 
an integration time of 50 ns. 

We note that an even better SNR can be obtained with 
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higher powers and a higher modulation index. However, 
as we shall see in Sec. IIVBI the presence of laser drift 
suggests that we choose the above values. The fact that 
our Pr'^+:Y2Si05 crystal has a diameter of 19 mm allows 
for a relatively high power without burning too deep holes 
causing drift problems. 

IV. LASER CHARACTERIZATION 

In this section we have two main objectives: to char- 
acterize the performance of our stabilized laser and to 
present some experimental verification of the theory con- 
cerning the drift calculations in Sec. Ill El 

The conventional method of characterizing laser per- 
formance is to build two identical lasers, setting up an in- 
terference experiment with these, and calculate the Allen 
deviation of the observed intensity beatings on different 
timescales. Since we only have a single stabilized laser 
we cannot apply this simple and useful method. Instead, 
we characterized the long-term drift by burning a spec- 
tral hole in an auxiliary crystal and measuring how the 
shift in laser frequency relative to this hole. For short 
timescales (up to the optical coherence time, T2, of the 
atoms) we use optical free induction decay (FID) to mea- 
sure the phase stability. We will not present experimental 
data on intermediate timescales, but the extremes above 
correspond exactly to the important cases for our appli- 
cations in quantum information processing. 

As already discussed in Sec. IIIFi there are three 
ground and three excited states in Pr^+:Y2Si05 , as 
shown in Fig. Hfe). The transition strengths between 
these are very different; the strongest and most con- 
tributing ones are ±i ±i (0.55), ±i ±| (0.38), 
±1 ^ ±i (0.40), ±1 ^ ±f (0.60), and ±| ^ ±f 
(0.93), while the remaining four are weak (0.07 or less). 
The numbers in parentheses are the relative strengths 
taken from flO] and we assume that these are also valid 
for branching ratios in the decay process. We see that the 
ions typically stay within the ±i, ±| space or in the ±| 
state, and only seldom change between these (the total 
crossing probability being 7%). If this crossing occurs, 
the RF pumping on the ground state hyperfine transi- 
tion ±1 ±1 at 17.31 MHz transition will counter- 
act it. Since the crossing is infrequent, the timescale 
for this RF transition Titmhz can be relatively slow (in 
order to maintain a certain hole depth). On the other 
hand, an ion resonant on, e.g., the ±^ ^ ±^ optical 
transition will decay to the ground ±| state with a high 
probability (40%), and the RF pumping on the hyperfine 
transition ±| ±i at 10.19 MHz must counteract this 
with a relatively short timescale Tiomhz (to maintain the 
same hole depth for this atomic species). In all experi- 
ments (apart from some of those illustrated in Fig. [TU]) 
we have set Titmhz = 10-Tiomhz- Hence, by combin- 
ing a small branching ratio for the optical decay taking 
an atom away from the resonant transition with a long 
timescale for the RF pumping bringing the atom back to 




FIG. 6: Experimental setup for characterizing the laser per- 
formance. Apart from the locking beam needed for the laser 
stabilization system we place an additional probing beam for 
characterizing the locking itself. AOM 1 is in double pass 
configuration and allows us to scan the laser beam frequency 
without any beam motion. AOM 2 allows us to shift back to 
the original stabilized laser frequency to characterize the holes 
when the laser is locked. An extra crystal in another cryostat 
is used for measuring the laser frequency drift on long time 
scales. 



the resonant transition (or vice versa), we ascertain that 
the parameter R is similar for the atomic species reso- 
nant on the five transitions mentioned above (according 
to the results in Tab. |T|). Since, also, the strengths of 
these five transitions are not that different, we expect to 
see experimental results not too different from the simple 
three-level model of Fig. |DJb). 



A. Experimental setup 

In addition to the laser stabilization system itself, the 
experimental setup used to characterize the laser perfor- 
mance is shown in Fig. [S] We have the possibility to 
send part of the stabilized beam through two acousto- 
optical modulators (AOM 1 and 2 in Fig. [SI both of 
which are A. A. Opto-Electronique modulators). The 
first one (AOM 1) has center a frequency of 200 MHz 
and a bandwidth of 100 MHz, and is placed in double- 
pass configuration. AOM 2 has a center frequency of 
350 MHz and a bandwidth of 200 MHz, and is placed 
in single-pass configuration. We may use the zeroth- or 
first-order diffracted beams from AOM 2 for experiments 
on the locking crystal or on an auxiliary crystal (5 mm 
thick 0.005% Pr^+iYaSiOs cooled to 2 K in an Oxford 
Spectromag cryostat). Both AOMs are driven by a Tek- 
tronix AWG520 1 GHz arbitrary waveform generator. 



1. Drift measurement setup 

In the first part of the long term drift experiments we 
repeatedly scanned AOM 1 in the frequency range 170 
MHz to 210 MHz, which in double-pass configuration 
becomes an 80 MHz scan from 340 MHz to 420 MHz. 
At the same time, AOM 2 is kept fixed at 350 MHz and 
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the minus-first-order diffracted beam is sent to tfie lock- 
ing cryostat. We thus scan the probing laser from —10 
MHz to 70 MHz, and in this way we are able to measure 
the absorption around the carrier hole at zero frequency 
and one of the side bands at 50 MHz. The beam is ex- 
panded to almost cover the entire crystal to probe all 
active ions and to obtain a weak intensity. In addition, 
the scan is fast (800 |xs), and accordingly the probing 
beam does not affect the hole structure. The scan is 
averaged over approximately 100 runs. After the lock- 
ing crystal the beam hits a detector (denoted "Det 2" in 
Fig. which is a Thorlabs PDB150A detector set to a 
bandwidth of 5 MHz. An example of this scanning signal 
is shown in Fig. [7] (green trace). Looking at the insets 
we clearly see ringing effects caused by the fast scan. It 
is possible to deconvolve this ringing back to the original 
absorption profile (black trace) by the method described 
in [45I, [4^. The scan speed, together with the detec- 
tion bandwidth, allows us do resolve structures as narrow 

as Sf = scan speed/bandwidth = ^'^^^^^^^ = 20 kHz. 
This method allows us to compare the center and side 
holes simultaneously, as shown in Fig. [lpra,b) . It should 



45l . |46| the absorp- 



be noted, that in the calculations in 
tion is assumed to be low {aoL ^1), and hence in our 
case corrections may be necessary to the reconstructed 
absorption. Since we already have a theory that is not 
exactly correct when approaching a^L 1 we shall not 
pay further attention to this fact. 

In the second part of the long-term drift experiments 
we used both the zeroth- and minus-first-order diffracted 
beams from AOM 2, and the scan width was narrower. 
The zeroth-order beam is used to first burn a spectral 
hole in the extra crystal (AOM 2 is turned off while do- 
ing this in order not to disturb the locking crystal) and 
the absorption from this hole is subsequently measured 
several times (with AOM 2 turned on). At the same 
time, the minus-first-order beam from AOM 2 measures 
the absorption profile for the center hole in the locking 
crystal. This setup enables us to measure the laser fre- 
quency drift and correlate the drift rate with the shape 
of the center hole. Examples of holes read out repeatedly 
for drift measurements are shown in Fig. [51 the detailed 
analysis of this will be discussed in Sec. IIVBI 



2. Phase stability measurement setup 

In order to measure the laser stability on short 
timescales we modified the setup shown in Fig. [5] slightly 
such that the zeroth-order diffracted beam from AOM 
2 (which is not used) is sent to the locking crystal with 
a beam diameter of roughly 2 mm. AOM 1 is operated 
around its 200 MHz center frequency and hence, with 
the frequency shifted around 400 MHz in double pass 
configuration, the probing beam will not interfere with 
the laser locking system. 

We used optical FID to measure the laser stability and 
to this end we programmed the pulse sequence shown in 
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FIG. 7: (Color online) a) The hole spectrum of the locking 
crystal, shown in black, observed after deconvolving the raw 
transmission, show in green. The peak labeled "1" and shown 
in inset b) at zero frequency corresponds to the carrier beam. 
This beam gives rise to several holes and anti-holes in a fre- 
quency range of ±36.9 MHz; the outermost one is labeled "2". 
This illustrates that a modulation frequency of cOm = 27r • 50 
MHz is desirable in order that the side hole, labeled "3" and 
shown in inset c), does not interfere with the hole spectrum 
of the carrier beam. The measured line-widths of the peaks 
shown in insets b) and c) are 60 kHz and 36 kHz for the carrier 
and sideband holes, respectively. The raw spectrum spanning 
-10 MHz to 70 MHz was collected in 800 ^s. 



Fig. M,c) and discussed in the figure caption. Referring 
to this figure, when pulse "2" is applied a coherence is set 
up in the atomic medium and when pulse "2" is turned 
off the atoms will keep radiating for a time limited by the 
optical coherence time, T2, (which in our case is around 
18 lis) and also by the inverse bandwidth of the actual co- 
herence. This decaying radiation "3" gives a fingerprint 
of the phase of the laser during pulse "2" . At the same 
time, we apply another pulse, "4" , shifted 45 MHz in fre- 
quency carrying its own phase. The beating of pulses "3" 
and "4" hence compares the present phase and the past 
phase, and this allows us to calculate the characteristics 
of the laser, as discussed in detail in Sec. IIVCI 

Fig.[5]illustrates a simple and useful method [l^l of ob- 
taining the amplitude and phase of an oscillating signal, 
in our case the FID heterodyne signal. The raw detector 
signal is Fourier transformed and filtered in a 40 MHz 
bandwidth around the positive 45 MHz component. The 
fact that the negative frequency components are removed 
leads (after a subsequent inverse Fourier transform) to a 
complex representation of the filtered FID signal in the 
time domain. The phase and amplitude can be read out 
directly from this complex signal. As an illustration, we 
show the filtered version of the FID signal in Fig. [H{a,b) 
together with a 45 MHz local oscillator signal derived 
from the 1 GHz arbitrary waveform generator. We see 
that the signals remain in phase over a period of 10 \is. 
This kind of measurement is repeated 100 times, which 
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FIG. 8: (Color online) Illustration of the phase stability mea- 
surements. In the main figure the lower (blue) trace shows the 
raw detector signal from the heterodyne FID detection. The 
upper trace (red) shows a filtered version of the lower trace. 
The lower traces (red) in insets a) and b) show this filtered 
signal in a 200 ns window around time and 10 (is, respec- 
tively. This can be compared to the local oscillator shown by 
the upper traces (green). We were able to obtain a good SNR 
for the phase difference of these signals. Inset c) shows the 
pulse sequence for the experiment: (1) is a 10 ms burn pulse 
scanning between 40 MHz and 50 MHz (relative to the AOM 
double-pass center). After waiting 100 ^.s, pulse (2) with a 
constant frequency of zero and duration 40 \is sets up a co- 
herence in the atoms which leads to the FID at (3). Finally, 
pulse (4) at a frequency of 45 MHz beats with the FID signal, 
leading to the detector signal shown in the main figure (lower, 
blue trace). 



allows us to calculate the statistics of the laser stability 
in a quantitative manner. 



B. Long-term drift 



Let us now turn to the experimental results concerning 
the laser drift and the hole shapes in the locking crystal. 
In Sec. Ill El we argued that if the spectral holes used 
for locking are too deep, the laser may be locked but 
drifting linearly in frequency, which in turn will cause 
the hole shape to be asymmetric. This is illustrated in 
Fig. fWa.b). From these measured center and side hole 
shapes, Q;R,(tt')L, we can calculate the imaginar y p arts, 
ai{uj)L, by using the Kramers-Kronig relations [48|. In 
our case, these relations take a slightly simpler form than 
usual since in Eq. ^ the only cj-dependence is in the 
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FIG. 9: (Color online) Example of drift measurements, a) 
Total read-out time is 3 seconds and we see a drift rate of 
roughly 160 kHz/s, while the direction changes occasionally, 
b) Total read out time is 20 s and the drift is 0.3 kHz/s over 
this time. 



denominator ^ + i{A — lu). It can be shown that: 
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(w - wo)2 -I- (52 



(61) 



Since cki(wo) is an integral of Q!r(cl;) times an odd func- 
tion in w — Wo I the value of Q!i(Ao) may be regarded as 
a convenient measure of the hole asymmetry, which al- 
lows us to quantitatively compare the center and side 
hole shapes. This is shown in Fig. llOf c) for a number of 
different settings of Tiomhz and Tnuuz- We see a clear 
proportionality, a'f""'\Ao) = 0.72 • af^'^-^Vo)- 

We expect aj'^'''™'^^ (Aq) and aj^"^"-' (Aq) to be equal, 
since from Eq. ([6]) the phase shifts, of the carrier (pc 
and sideband 0s are proportional to aj'^'*"'"''' (Aq) and 
^(sido) j^^^-j^ respectively, and with a closed laser stabiliza- 
tion feedback loop we must have zero error signal with 
0c = 0s- The reason for the slope not being unity is 
unknown. We have thus shown that the laser may drift 
linearly, while the feedback loop is still locked. We base 
this on the direct observations of the drift, as exemplified 
in Fig. [Qlja), on the proportionality in Fig.fTOlfc). and on 
the fact that the electronic error signal is small. 

We may also wish to correlate the measured hole 
asymmetry, q;i(Ao), with measured drift rate, (3, since 
these should be related by Eq. ((57)) . at least for small 
drift rates when the first-order theory is valid. This 
equation assumes a three-level model, and for our real 
Fr^+ : Y2 SiOs ions we observe that the±i,±|^±i,±| 
transitions with typical branching ratios around 0.5 and 
hole lifetimes, given mostly by Tiqmhzj are reasonably 
representative. Hence, we use T^g = Tg^ = Tiqmhz and 
6or = 0.5. Also, from the experimentally measured hole 
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FIG. 10: (Color online) In the two upper graphs the measured 
absorption, q:r(cj)L, is plotted for comparison under different 
conditions for a number of center holes (a) and side holes (b). 
With increasing hole depth, the asymmetry of both increases. 
The imaginary part, qi(Ao), of the absorption length at the 
hole center is a quantitative measure of this asymmetry, and 
can be calculated from aR(tj) using the Kramers-Kronig re- 
lations l|61|) . In (c) we see, for several measurements under 
different conditions, a clear linear relationship between this 
asymmetry for the center and side holes. The straight line is 
a fit through the origin with a slope of 0.72, theoretically we 
expect a slope of unity. 
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we can estimate 

Fhoic and dhole from Eq. (|22|) . Inserting these parame- 
ters together with the experimentally determined drift 
rate, /3, into Eq. ([57|) we can calculate, to first order, the 
value of ai(Ao) for the center hole. This is plotted on the 
abscissa in Fig.ITT] On the ordinate Q!i(Ao) is plotted, in- 
ferred by the Kramers-Kronig relations in Eq. (|6ip from 
the measured aR(ijj). The result shows a proportional- 
ity between these two and hence we have demonstrated 
a clear correlation between measured hole shape asym- 
metry and measured drift rate as suggested by Eq. ([57)1 . 
The fact that the constant of proportionality is not unity 
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FIG. 11: (Color online) Comparison of the measured and 
calculated imaginary center hole absorption depth, ai(Ao)L, 
for different settings of the RF pumping times Titmhz ~ 
10 ■ TioMHz. On the y-axis ai{Ao)L, as inferred from the mea- 
sured aB.{u!) via the Kramers-Kronig relations in Eq. (I6ip 
is plotted. On the x-axis ai(Ao) calculated from Eq. (|57|) 
with the experimentally inferred values of aoL, Fhoie.c, dhoic.c, 
Trg = Tgr = TioMHz, and /3, while Ti = 150 \ls and 6er = 0.5 
is plotted. The blue circles and red squares correspond to the 
data shown in Fig. I12f a') and Fig. I12f c'). respectively. These 
two cases differ by 1 GHz in the position on the inhomoge- 
neous profile, giving background absorptions, oqL, of 0.44 and 
0.66 for the blue circles and red squares, respectively. 



is not alarming since we used a simplified model with 
only three atomic levels. Also, when the laser is drifting 
the first-order theory should not be sufficient. The figure 
shows results from two different experimental runs at two 
different positions on the inhomogeneous profile. These 
have a different slope, which is not understood. From 
the results in Fig.[TT]we have demonstrated that the pre- 
dictions of Eq. ([57)1 show the correct order of magnitude 
when compared to experiment. It should be borne in 
mind that a full understanding of the data requires more 
than our first-order calculations. 

Let us now turn to the measured drift rate, /3, for var- 
ious parameter settings. The results of these measure- 
ments are shown in Fig. 1121 The light intensity is kept 
constant with the saturation parameter sq ~ 0.09. The 
modulation index, m, has the values 0.14, 0.20, 0.28, 0.40, 
and 0.56, and each color in the figure corresponds to one 
of these values. The hole shapes are controlled by em- 
ploying TioMHz = jq ■ TirMHz at the six values 2 ms, 4 
ms, 8 ms, 20 ms, 40 ms, and 80 ms. In the figure the 
six corresponding data points are plotted for each color 
from left to right since the hole depth <i|™i°''c^ increases 
with TioMHz. 

The results in Fig. [T^] are divided into two parts. In 
Fig. [T2Ua) we see that the drift rate is lowest for hole 
depths around 0.5, corresponding to Tiqmhz = 8 ms. For 
shorter timescales (towards the left) the drift rate tends 
to increase and this could be caused, for example, by off- 
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FIG. 12: (Color online) Measured drift rates versus measured 
hole depth. The hole shapes are changed by varying Tiomhz 
and TiTMHz- Red circles, m = 0.56; green crosses, m = 0.40; 
blue squares, m = 0.28; cyan triangles up, m = 0.20, purple 
triangles down, m = 0.14. The vertical dashed lines indicate 
the values of the corresponding thresholds shown in Fig. [S] 
The two graphs differ in that the data shown in (a) were 
collected before (c) and between the measurements the laser 
frequency was moved 1 GHz on the inhomogeneous profile, 
changing the background absorption, aoL, from 0.44 to 0.66. 
In addition, as shown in inset (b), the local oscillator phase 
(relative to the value used in (a)) to demodulate the error 
signal, was adjusted in order to minimize the drift rate for 
shallow holes. 



that the drift rate for a shallow hole depends on the po- 
sition on the inhomogeneous profile, which was then ad- 
justed by 1 GHz. In addition, as shown in the figure inset 
(b), adjusting the local oscillator phase slightly also has 
an impact on the drift rate, which we fine tuned to a low 
value. Fig. [TWc) shows low drift rates in the left part of 
the figure (the lowest measured being below 0.5 kHz/s). 
However, the higher drift rates in the right-hand part of 
the figure remain almost unchanged. A slight change in 
error signal offset cannot change the fact that zero drift 
is an unstable solution if the criterion in Eq. ([55)1 is not 
met. The increase in drift rate on the right-hand side of 
Fig. I12f c) is now consistent with the vertical lines repre- 
senting the threshold. The highest drift rate is still below 
25 kHz/s which is fairly good. 

A very rough order of magnitude estimate of the drift 
rate can be obtained by noting the fact that the calcu- 
lations in Sec. Ill El rely on the series expansion in the 
dimensionless parameter, ^ = p . To calculate the 
threshold condition for (3 = bemg a stable solution it 
is sufficient to assume that ^ <C 1. However, to calcu- 
late the actual drift rate when /3 = is unstable requires 
higher order theory, where ^ cannot be small. Thus, if 
we then take for our right-most data points the typical 
values Trg = 80 ms, Thoic/ZTr « 100 kHz, and /3/27r = 25 
kHz/s, we obtain ^ « 0.02, which is actually quite low 
compared to unity. This is an indication that the fraction 



of ions resonant on the ± | 



± I transition with ten- fold 



longer hole lifetime, Trg = TnMUz, and ten- fold higher 
^ plays a role in limiting the measured drift rate to a 
low value. Without having completely searched the very 
large multi-dimensional parameter space, our experience 
tells us that Tiomhz — jq ■ Titmhz is a good choice. We 
have observed examples of much higher drift rates than 
25 kHz/s, see e.g. Fig. IDJa) where, in addition to a high 
drift rate of 160 kHz/s, we also saw that the direction 
changed occasionally. 



set errors in the electronics, or by the small signal from 
the inhomogeneous background if the phase of the local 
oscillator is set incorrectly. Note, that if the error sig- 
nal has an offset of 10"'^ times the full-scale value, we 
would expect the laser to be displaced by 10~'^rhoio in 
frequency. With a hole lifetime of T^g we can then es- 
timate the drift as /3 « lO^'^Fhoic/Trg. For our shortest 
timescales T^g « 2 ms, we have Fhoic/27r w 40 kHz lead- 
ing to an estimate of P/2tt « 20 kHz/s. This is not far 
from the experimental values, which we find to be typi- 
cally a little less than 20 kHz/s and we thus concluded 
that our relative offset errors are smaller than 10"'^ for 
the data shown in Fig. [T^ a). For longer timescales (to- 
wards the right) we can also see an increase in drift rate, 
which can be explained by the fact that we are approach- 
ing the maximum hole depth, discussed around Eq. (j59p 
and Fig. [31 For each color in Fig. [T^] this threshold is 
shown as a vertical dashed line. 

The above conclusions become more apparent when we 
consider Fig.fTWc). Before collecting these data we noted 



C. Short-term phase stability 

Fig- flWa.b) shows the phase evolution of the filtered 
FID signal for the case when the laser stabilization sys- 
tem is on (black) and off (red). In the latter case, the dye 
laser is stabilized to a reference cavity by its commercial 
feedback system. The clear effect of our locking system 
can be seen, and below we analyze this phase evolution 
quantitatively. 

We note from Fig. [8] that our method based on de- 
tection of the FID signal works in practice for times 
up to around 15 [is. The SNR when determining the 
phase will decrease with increasing time. Also, referring 
to Fig. Efc), the variation in the phase of the laser is 
present in both the FID signal "3" and in the heterodyn- 
ing pulse "4" . In order to derive the actual phase stabil- 
ity of the laser we repeated the measurements shown in 
black in Fig.fTWa') 100 times, divided the measured phase 
into time bins of length 500 ns, and performed statistical 
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FIG. 13: (Color online) (a) A typical single-shot example 
of the measured FID phase with our stabilization system on 
(black) and off (red), (b) An enlargement of the first ii.s of 
the plot shown in (a), (c) Phase standard deviation versus 
time, collected in time bins of 500 ns. Red diamonds repre- 
sent the experimental uncertainty in the measurements, given 
by Eq. (|62|l . Green circles show the shot-to-shot variations, 
std(a;i). The black circles describe the variation of the phase 
with respect to the phase of the first time bin, std(a;i — xi). 
The green and black squares show the results when the mea- 
surement noise has been subtracted. The blue line shows the 
absolute value of the mean phase compared to the first time 
bin. The black line is a linear fit to the black squares, giving 
(5(^[deg] = 2.0 + 2.0-10^t[s]. 



operations on these measured phases. Note, that when 
discussing the phase of the FID signal we always mean 
the phase difference between the reference signal and the 
FID signal (shown in green and red in Fig. [5I^a,b), re- 
spectively). 

In the following, Xij denotes the measured phase in 
time bin i for repetition j = 1, . . . , 100. For each time 

'i — n 

In Fig. [TST c) abs(a;i) is shown by the blue line and 
std(a;i) = A/Var(xi) by the green circles. We see that 



bin i we may calculate the mean value Xi — Xij and 

variance Var{xi) — ~ ^^j)^, where n = 100. 



the average phase deviates less than 1 degree over time, 
which simply means that there is no phase drift in one 
particular direction. The std(xi) values describe the re- 
peatability of the experiment, which is within roughly 
5 degrees and increases with time. To characterize the 
measurement noise we now assume a model where the 
measured values can be written x 



^,3 



riij, where 



is the actual phase of the FID signal and T^i^j IS a 
stochastic variable describing the noise in the j'th time 
bin and j'th measurement. Since technical laser noise 
in general is very limited at the heterodyne beating fre- 



quency 45 MHz, we assume the noise level to be shot 
noise or electronic noise in the detector. This is broad- 
band (white) noise and we may assume all the individual 
are independent. The mean value is assumed to be 



equal to zero, and the standard deviation std(rti) will de- 
pend on i due to the decaying FID signal. In order to 
estimate the standard deviation of n^j we calculate the 
variance of adjacent time bins for our experimental data: 



Yai{xi - Xi-i) = Var((/)j - (t>i^i) + Var(ni 
~ const -I- 2 • Var(ni). 



l) 

(62) 



In the second step we use the fact that Var(0i — 0i_i) 
must be time independent since the laser phase is in a 
steady-state condition. We also use the approximation 
Var(ni — nj_i) — Var(nj) + Var(?T.i^i) sa 2Var(nj), since 
we assume the detection noise to vary slowly over time. 
Plotting the variance \a,v{xi — Xi-i) versus time (i.e. plot- 
ting Eq. ([62)1 as a function of i) will give the measure- 
ment noise, Var(ni), and the constant, Var(^i — 
In Fig. [TSFc) std(ni) = yj\ai{ni) is shown by red dia- 
monds, and is seen to increase from zero to less than 4 
degrees at long times. 

This noise level can now be used to correct the shot-to- 
shot variance, shown by the green circles, since Var(a;i) — 
Var(0i) -|- Var(ni). Subtracting the measurement noise, 
the green squares show the standard deviation, std((/ii), 
of the FID signal itself. 

When characterizing laser performance it is more in- 
teresting to measure the phase evolution over time, 
i.e. how much the phase of the i'th time bin, ^j, devi- 
ates from the first one, 0i. To this end we first calculate 
Var(a;i — xi) = Var((/)i — 0i) -I- Var(ni — ni) from the ex- 
perimental data. The square root of this is std(a;i — Xi), 
which is plotted in Fig. fTST c) as black circles. Since 
Var(ni — ni) = Var(ni) -I- Var(ni) — Var(ni), we can sub- 
tract the measurement noise again to obtain std(0i — (/)i), 
which is plotted as black squares. We see that on short 
timescales the phase error is around 2 degrees, and it in- 
creases linearly with time. The black line in the figure 
is a linear fit to the black squares, and it can be used to 
convert the phase stability into a characteristic frequency 
stability as a function of time. Since the relation between 
phase and frequency is </'[deg] = 360 • / [Hz] - t + (p^, the 
standard deviations of these will be connected by: 



std(/) = 



std((/)) _ 2.0 -f 2.0-104 



360 -i 



360 



(63) 



In the second step, we inserted the experimental linear 
fit shown as the black line in Fig. fTSTc). For large t 
the characteristic frequency standard deviation, std(/), 
approaches 0.6 kHz, but we are not allowed to extend t 
further than the measured 14 (a.s. Note, for t = 10 ps, 
std(/) = 1.1 kHz. 

The difference between the shot-to-shot standard de- 
viation, plotted in green in Fig. I13r c). and the time sep- 
aration standard deviation, shown in black in Fig. I13f c). 
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is approximately 3 degrees. The extra noise in the shot- 
to-shot data arises from the fact that the laser phase also 
fluctuates during the setup of the coherence in the atoms 
(pulse "2" in Fig. [5]). The practical phase memory time is 
approximately 10 ij.s, corresponding to the decay time of 
the FID signal in Fig. [H Averaged over 10 \JlS, the black 
data points of Fig. [T3] actually suggest approximately 3 
degrees fluctuation in the laser phase. Thus the data 
shown in green and black in Fig. [13] are consistent with 
each other. 

The experimental data shown in Fig. I13r c) were ob- 
tained with TioMHz = jq ■ T'lTMHz = 4 ms, modulation 
index m — 0.20, and saturation parameter sq « 0.09. 
This setting is also represented in Fig. [T2Fc) giving a 
drift of (0.34±0.76) kHz/s. Hence, we have shown that a 
slow drift and low phase variations can be obtained at the 
same time. Of the settings investigated the best general 
phase stability performance is found for the above RF- 
eraser time, which gives rise to a measured hole depth, 
4Tc'!c^ around 0.5 to 0.6 in Fig. [H This hole depth 
corresponds to x = « 2, which is the best choice 

according to the /-function ([38]) shown in Fig. [1] At 
other settings we still flnd asymptotic values of Eq. ((63|) 
less than 2.5 kHz, which is only four times worse than 
the value of 0.6 kHz in the example above. 

We have also investigated a few different positions on 
the inhomogeneous profile. These measurements indicate 
that a transmission above 30% is a good choice. Varying 
in the local oscillator phase between the values shown 
in Fig. [T2jb) does not change the performance on short 
timescales. 

The performance of the non-stabilized laser, i.e. when 
the dye laser is operated only with its built-in stabiliza- 
tion unit, is exemplified in Fig. [T3ja,b). We see that 
in this case it takes roughly 13 \JlS to change the phase 
by 360°, which can be translated into a characteristic 
frequency of 75 kHz. Repeating this measurement 100 
times shows that these values are typical. In our experi- 
ence, cooling the laser dye affects this frequency consid- 
erably (we typically use temperatures in the range 8°C 
to 12°C). 

The above observation does not mean that the unsta- 
bilized laser line-width is 75 kHz, it only means that on 
short timescales of the order of lOus, the frequency de- 
viations are around 75 kHz. In [l9| a typical dye laser 
line-width of 450 kHz is reported, but it was also in- 
dicated that the frequency deviations had a bandwidth 
B <C 190 kHz, i.e. the deviations occurred on a timescale 
much longer than 5 |J.s. Our observations are consistent 
with this behavior. 

As discussed in Sec. IIIIl the group delay of our sta- 
bilization feedback loop is approximately 100 ns. To 
understand how this limits the phase stability of our 
laser we assumed that the dye laser introduces a fre- 
quency error of 75 kHz. It takes 100 ns for the sta- 
bilization loop to discover this and start counteracting 
it, and after this time the phase has evolved roughly 
360° • 75 kHz • 100 ns = 2.7°. This hand- waving argu- 



ment is consistent with a standard deviation of the order 
of 2 degrees, as shown by the data (black squares) in 
Fig. fTST c) as the limiting case for short timescales. An- 
other way of saying this is that in Fig. [T^ b) the typical 
stabilized phase deviation (shown in black) corresponds 
to the actual deviation of the non-stabilized laser (shown 
in red) after a time of the order of the delay time in the 
feedback loop. 



D. Overall laser performance 

In summary, we can say that we achieved 1 kHz line- 
widths on the timescale of 10 |j.s, and drift stabilities 
also around 1 kHz on 1 s timescales. We were not 
able to thoroughly examine the stability on intermedi- 
ate timescales. If we had had two identical lasers with 
which to calculate the Allen deviation from beating ex- 
periments, our reported line widths should be multiplied 
by \/2 for short timescales (with random shifts) and by 
2 for long timescales (with linear drift). Then we can 
compare our laser performance with others reported in 
the literature. 

Laser stabilization using a spectral hole has been re- 
ported for semiconductor lasers in several publications 
[H, m, m, m, m, [23, Ull. in these cases the reported 
Allen deviation on a lOus timescale is typically 10 kHz, 
and around 3 kHz in [23|. For long timescales a linear 
drift rate of 10 kHz/s has been observed in pure tran- 
sient systems [23] • Stabilities better than our reported 1 
kHz/s can be obtained with permanent spectral holes (25 
kHz/min) ^2^, (« 1 kHz/min) \2§\ or by incorporating 
the permanent reference consisting of the inhomogeneous 
|)rofile (7 kHz/min) 23], (1.4 kHz/min, no hole burning) 

The most important steps in obtaining our short-term 
phase stability of down to 2 degrees, consist of (1) con- 
structing fast analog electronics to obtain short delays 
and low noise, (2) using a large sample crystal diameter 
(19 mm) in order to be able to use a high light power 
without too high an intensity, and (3) cooling the laser 
dye in order to improve the passive stability. We should 
also mention that the optical table is placed on legs with 
pneumatic vibration dampers, that an improved dye jet 
nozzle (Radiant Dyes RDSN 02) and an improved nozzle 
holder and pump mirror adjustment unit (Radiant Dyes 
RDU 10) have been installed, and that the entire laser 
system is located inside a clean room (from Terra Uni- 
versal). In order to obtain the low frequency drift of 1 
kHz/s the important steps are: (1) being aware of the 
drift problem discussed in Sees. Ill El and IIVBI and (2) 
carefully constructing the electronics so as to ensure low 
offsets. 
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1. Impacts on quantum gate experiments 

As we discussed in the introduction in Sec.[Tl our main 
reason for frequency stabilizing the laser is to enable 
quantum computation experiments. One experimental 
approach utilizes spectral hole burning to isolate strongly 
interacting qubits, as described in [To, [HI. Such qubit 
structures will typically have a spectral width of ~ 200 
kHz, and can be prepared by several hundred cycles of 
optical pumping. The entire preparation process may 
take hundreds of milliseconds, and it is important that 
the laser drift is much less than 200 kHz during this time. 
With typical drift rates of 1 kHz/s we have fulfilled this 
requirement. 

Let us now consider the effects of laser phase fluc- 
tuations on the fidelity of qubit operations. We use 
two ground state hyperfine levels to represent a qubit, 
e.g. |0) = jij) and |1) = |±f). Single qubit operations 
can be performed by interaction with an excited state 
|e) [1] where 7r-pulses are used to transfer the popula- 
tion from |0) or |1) to |e) and back again. Such basic 
operations are also building blocks for more complicated 
gates, e.g. the CNOT gate [l|, Here we will estimate 
the impact of a phase error in the laser on such a simple 
operation. 

Consider the following scenario. The qubit is in an 
initial state |V'o) = ct \0)+f]\l), then a perfect laser drives 
a TT-pulse on the |1) — > |e) transition such that \iI!q) 
a\0) + (3\e). The laser phase is now assumed to change 
by an amount (j), and we finally drive the excited state 
back to the final state ^ a\0) + /Je'*^ We have 
thus gathered all the phase fluctuations in a single step. 
If there is no phase change (0 = 0) we arrive back at the 
initial state \tpo) , and the deviation from this can be used 
to estimate the fidelity of the operation given (j). 

However, with an unknown random phase change, 0, 
the correct description of the final state is given by the 
density matrix: 

Pf= / |M^f|/(0)d0=|ap|O)(O|-f |/3ni)(l| 

^ r (64) 

+ J (a/3*e-'*|G)(l| + a*/3e**|l)(0|)/(</))d<^, 

where /(</>) is the probability distribution of the phase 
change, (p, where J f{(f))d(f> = 1. Given an initial state 
{(po) the fidelity of the operation is the overlap: 

F = (Vol PI l-Z-o) = \a\^ + |/?|' + 2\a\^\f3\^Jf{cf,) cos0d0 
= 1 - 2|an/3pe with e = 1 ~ J f{(j)) cos (l)d<j). 

(65) 

The parameter e is zero when there is no phase change 
if {4') = H4')) and varies to second order in the phase 
deviation. To see this clearly we assume a simple top- 
hat distribution function, f{(j))= if —\/3(f>o < (p < 



VS^o and zero elsewhere. This distribution function cor- 
responds to a root mean square phase (pQ. Inserting this 
in Eq. ^ gives e ^ \(Pl and F ^ I - \a\'^\(3\'^ (pl when 
00 < 1- 

The phase fiuctuation plays no role if the initial state 
is either of the basis states, |0) or |1). The impact is 
maximal if the initial state is an even superposition with 
lap = = ^. In the latter case, for a phase error with 
a standard deviation as high as 10 degrees, the fidelity is 
98.5%. 

More detailed analysis is naturally required to accu- 
rately calculate the fidelity loss of entire quantum gate 
experiments given our laser stability, which is outside the 
scope of this paper. However, with the above investiga- 
tions, we are confident that the laser stability is sufficient 
for our intended quantum gate experiments. 



V. OUTLOOK 

As discussed in Sec. IIVD 11 we have a laser in a suffi- 
ciently stable condition for quantum information exper- 
iments, which was the main reason for the entire stabi- 
lization project. Compared to previous laser stabilization 
studies on pure transient hole-burning systems, we have 
achieved better stability, and the understanding of laser 
drift was very important in this process. However, some 
compromises had to be made since the absence of linear 
drift does not necessarily coincide with the maximization 
of error signals. The presence of both the center and side 
holes is responsible for this effect, and a number of mea- 
sures can be taken. (1) An interferometric setup can be 
considered where only the beam not passing the atoms is 
modulated in the Pound-Drever-Hall scheme. This, how- 
ever, requires a very stable or an otherwise stabilized in- 
terferometer, since the error signal from the spectral hole 
locking can not conveniently distinguish whether phase 
changes occur relative to the atoms or because of insta- 
bilities in the interferometer. (2) Moving the modulation 
side bands outside the inhomogeneous profile. This is 
inconvenient in our case with Pr^+:Y2Si05 since the in- 
homogeneous broadening is around 5 GHz. However, the 
system under consideration in ^40] is interesting in this 
respect. Note, that the drift problem decreases when 
'^m/rinh approaches unity, so a combination of the inho- 
mogeneous profile for permanent stability together with 
hole burning effects for short-term stability may be pos- 
sible without suffering too strong constraints from the 
drift criterion when lj^ « Tinh- 

In any case, our work has also shown that offset lev- 
els and systematic effects must be carefully controlled to 
achieve low drift rates in pure transient systems. For this 
reason, the incorporation of permanent effects seems to 
be the most convenient and competitive solution for ob- 
taining a slow long-term drift. It is our hope that the 
work presented here will help to further improve the al- 
ready advanced field of laser stabilization. 
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VI. CONCLUSIONS 

We have for the first time stabiHzed a dye laser to a 
spectral hole in Pr"^+:Y2Si05 , obtaining 1 kHz frequency 
stability on a 10 |j.s timescale, and a long-term drift rate 
below 1 kHz/s. The use of RF fields to control the hyper- 
finc level populations allowed us to optimize the spectral 
hole parameters for best laser performance. The stabil- 
ities obtained were sufficient for high-fidelity quantum 
information experiments. 

We have contributed to the theoretical understanding 
of laser stabilization using spectral holes to an extent 
that we hope will enable other scientists to further im- 
prove existing technology. Although a few experimental 
observations reported in this paper are not fully under- 
stood in a quantitative manner, we have provided strong 
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